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ABSTRACT 

The  classical  Stefan  theory  for  the  melting  of  a  solid  or  the  freezing  of 
a  liquid  is  too  simplistic  to  describe  phenomena  such  as  supercooling,  in 
which  a  liquid  supports  temperatures  below  its  freezing  point,  or 
superheating,  the  analog  for  solids,  or  dendritic  growth,  in  which  simple 
shapes  evolve  to  complicated  tree-like  structures.  In  this  paper  we  develop  a 
general  theory  for  two-phase  phenomena  of  this  type  starting  from  general 
thermodynamical  laws  which  are  appropriate  to  a  continuum  and  which  include 
contributions  of  energy  and  entropy  for  the  interface  between  phases.  We  show 
that  the  interfacial  temperature  is  generally  not  equal  to  the  melting 
temperature,  but  depends  in  a  prescribed  manner  on  the  curvature  of  the 
interface.  We  describe  appropriate  initial-boundary-value  problems,  at 
various  levels  of  approximation,  and  deduce  corresponding  Liapunov  functions. 
We  also  present  a  general  theory  of  equilibria,  and  show  that  -  for  an 
unbounded  domain  in  isotherma1  equilibrium  -  there  are  no  stable  states  in 
which  the  bounded  phase  is  solid  and  the  unbounded  phase  supercooled  liquid. 

We  show  further  that  corresponding  minimizing  sequences  for  the  free-energy 
are  consistent  with  the  formation  of  dendrites. 
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SIGNIFICANCE  AND  EXPLANATION 


— /  The  classical  Stefan  theory  for  the  melting  of  a  solid  or  the  freezing  of 
a  liquid  is  too  simplistic  to  describe  phenomena  such  as  supercooling,  in 
which  a  liquid  supports  temperatures  below  its  freezing  point,  or 
superheating,  the  analog  for  solids,  or  dendritic  growth,  in  which  simple 
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shapes  evolve  to  complicated  tree-like  structures.  In  this  paper  we  develops, a 

ti 

general  theory  for  two-phase  phenomena  of  this  type.  We-  develop^ partial 
differential  equations  satisfied  in  the  phase  regions  and  free-boundary 
conditions  satisfied  on  the  interface  between  phases,  and  give^arguments  which 


indicate  that  the  resulting  boundary-value  problems  predict  the  formation  of 
dendrites .  j  ^ 
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Morton  E.  Gurtin 

1.  Introduction. 

A  classical  problem  of  mathematical  physics  is  the  Stefan  problem  for 
the  melting  of  a  solid  or  the  freezing  of  a  liquid.  The  underlying  theory, 
however,  is  too  simplistic  to  describe  phenomena  such  as  supercooling,  in 
which  a  liquid  supports  temperatures  below  its  freezing  point,  or 
superheating,  the  analogous  phenomenon  for  solids,  or  dendrite  formation 
(cf.  Figure  1 ),  in  which  simple  shapes,  such  as  spheres,  evolve  to 
complicated  tree-like  structures.*  The  past  two  decades  have  seen  the 
development  of  more  general  theories2  for  phenomena  of  this  type,  a 
critical  ingredient  being  a  free-boundary  condition  at  the  solid-liquid 
interface  I  *  I(t)  in  which  the  temperature  depends  on  the  curvature  of  1. 

In  these  theories  questions  arise  as  to  what  are  the  interface  conditions;3 
in  fact,  it  is  not  clear  which  of  the  interface  conditions  are  constitutive 
assumptions  and  which  follow  directly  from  the  underlying  balance  laws. 

It  is  the  purpose  of  this  paper  to  develop  a  theoretical  framework  for 
theories  of  this  type  starting  from  general  thermodynamical  lav/s  which 
are  appropriate  to  a  continuum  and  which  include  interfacial  contributions 
for  both  energy  and  entropy.  We  do  not,  however,  seek  the  broadest  possible  theory, 
but  rather,  to  focus  our  attention  on  the  thermodynamics  and  to  keep  the  under¬ 
lying  relations  transparent,  we  limit  our  discussion  to  rigid  heat  conductors.4 
to  rigid  heat  conductors.4 

Our  chief  assumptions  -  apart  from  general  equations  of  state  for  the 
bulk  and  interfacial  quantities  -  are  that  the  interface  I  produce  no 
entropy  and  that  the  temperature  be  continuous  across  I.  Among  our  main 
results  are  the  interface  conditions 

(q]*m  *  vlE]  -  vfee  -  eA  on  I, 

T  *  ([E]  -  fee)/((S]  -  fes)  on  I,  (1.1) 

vm*n*0  on  dl, 

in  which  T  is  the  temperature;  IE],  [5],  and  Iq]  are  the  jumps5  in 
energy,  entropy,  and  heat  flux  across  the  interface;  e  and  s  are  the 
Interfacial  values  of  energy  and  entropy;  6*  v,  and  m,  respectively,  are, 
for  the  interface,  the  sum  of  principal  curvatures,  the  normal  velocity,  and 
a  unit  normal  vector  (with  sign  convention  explained  in  the  text);  eA  is 

*Cf.  Chalmers  1 15641  and  Delves  1 1974J  for  discussions  of  these  phenomena. 

2f1ultins  and  Sekerka  (1963.1964],  Voronkov  (19651.  See  also  the  review  articles  by 
Sekerka  ( 1968,1973.1984],  Chernov  (1972],  Delves  (1974],  and  Longer  (19SO ]. 

3Cf.  Rogers  1 1933]  for  a  discussion  of  soma  of  the  inconsistencies  in  the  literature. 

future  paper  will  discuss  the  effect  of  varying  concentration. 

50ur  convention  for  jumps  and  for  the  latent  heat  L  is  ‘phase  2  minus  phase  1*.  with  phases  labeled 
bo  that  L  i  0.  Thus  for  a  solid-liquid  system  phase  2  would  denote  the  liquid,  and  for  that  reason  we 
will  often  refer  to  phases  1  and  2.  respectively,  as  the  solid  and  liquid  phases.  _ 

Sponsored  by  the  United  States  Army  under  Contract  Nos.  DAAG29-80-C-0041  and 
DAAG29-82-K-0002  and  by  the  National  Science  Foundation  under  Grant  No. 
DMS-8404116. 
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the  time  derivative  of  e  following  the  interface;  n  is  the  outward  unit 
normal  on  the  boundary  of  the  region  B  occupied  by  the  body. 

The  first  of  (1.1)  is  essentially  the  first  law  of  thermodynamics  at  the 
interface.  The  second  -  derived  within  the  fully  dynamical  theory  -  is  a 
condition  of  local  equilibrium  expressing  balance  of  free-energy  across 
the  interface.  The  third  is  a  contact  condition  for  that  portion  of  the 
interface1  which  intersects  the  boundary  of  B;  it  assserts  that  -  where 
the  interface  meets  dB  -  it  is  orthogonal  to  dB  or  stationary. 

We  discuss  two  types  of  boundary  conditions:  an  isolated  boundary  on 
which  qn  =  0;  an  isothermal  boundary  on  which  T  is  constant.  We  show 
that,  for  either  of  these  boundary  conditions, 

interfacial  area  is  uniformly  bounded  in  time,  ( 1 .2) 

at  least  when  B  is  bounded 

We  also  discuss  the  equilibrium  theory  under  Isothermal  boundary 
conditions,  and  define  stable  states  as  minimizers  of  a  global  free-energy. 
We  consider  a  material  for  which  the  bulk  free-energies  cross  at  a  single 
temperature  TM,  and  show  -  for  bounded  B  -  that  stable  states  are 

always  single  phase,2  the  stable  phase  being  the  phase  with  lower 
free-energy.  Tn  therefore  represents  the  temperature  at  which  a  change 

in  stable  phase  occurs,  and,  for  that  reason,  we  refer  to  TM  as  the 
transition  temperature. 

The  question  of  stability  for  unbounded  B  is  far  more  interesting. 

Here  our  results,  expressed  in  terms  of  a  solid-liquid  system  in 
isothermal  equilibrium,  assert  that:3 

1.  There  are  no  stable  states  in  which  the  bounded  phase  is  solid  and  the 
unbounded  phase  supercooled  liquid. 

2.  Under  the  conditions  of  ( 1 ),  minimizing  sequences  of  the  free-energy 
are  consistent  with  interfacial  instabilities  such  as  the  formation  of 
complicated  arrays  of  thin  spikes,  behavior  indicative  of  dendritic  growth. 

Next,  we  Introduce  a  quasi-static  model  for  situations  in  which  the 
interface  moves  slowly  compared  with  the  time  scale  for  heat  conduction. 
The  chief  constitutive  hypothesis  underlying  this  model  is  that  -  in  each 
of  the  phases  -  both  the  bulk  energy  and  the  bulk  entropy  are  constant.  We 
also  assume  that  the  conductivities  k{,  the  interfacial  energy  e,  and  the 

1  Here  it  is  tacit  that  JlcdB. 

2Here  it  is  important  to  emphasize  that  the  boundary  is  held  at  constant  temperature:  two-phase 
solutions  are  possible  when,  for  example,  the  body  is  isolated  and  the  total  energy  constrained  (cf. 
Section  3.3). 

3ln  their  pioneering  paper  of  ( 1963},  Mullins  and  Sekerka,  working  within  the  dynamic*!  theory 
described  by  ( 1.13),  established  the  instability  of  the  interface  for  infinitesimal  perturbaUons  of  a 
sphere  solidifying  in  a  supercooled  melt.  The  assertions  (1)  and  (2)  are  analogs,  within  the  equilibrium 
theory,  of  the  Mullins-Sekerka  instability. 


interfacial  entropy  s  are  constant.  Let  Bj(t)  denote  the  subregion  of  B 
occupied  by  phase  i  (i  =  1,2),  and  let 


Then  these  assumptions  lead  to  the  system1 


au  *  0,  q  =  -k^gradu 

in  Bj, 

u  -  -  hfe/(  1  -afc),  [q]  m  =  (L  -  fee)v 

on  I. 

(1.3) 

vmn  =  0 

on  dl, 

where 

h  =  T„f(TM)/L,  a-TMs/L 

(1.4) 

with  f(0  the  interfacial  free-energy  and  L  the  latent  heat.  We  are  able 
to  establish  global  growth-conditions  for  (1.3)  under  the  two  types  of 
boundary  conditions  discussed  previously.  In  particular,  letting 

U  =  T„f(T0)/L,  O  •  e/L, 


we  show  that: 

(1)  for  an  isolated  boundary, 

voKBj)*  +  Barea(I)*  *  0,  area(I)*  i  0;  (1.5) 

(ii)  for  an  isothermal  boundary, 

u0vol(B2)*  ♦  parea(I)*  i  0.  (1.6) 

A  detailed  analysis  of  (1.5)  and  (1.6)  Is  beyond  the  scope  of  this  paper. 
We  do,  however,  give  a  lengthy  dlscussslon  and  indicate  several  results 
assuming  the  stability  of  solutions  as  t-*-«>.  In  particular,  we  give 
arguments  in  support  of  the  interfacial  instabilities  described  in  (1)  and 

(2)  above.  The  result  (1.5)  seems  also  to  indicate  an  instability 
characterized  by  a  solid  phase  whose  volume  tends  to  zero,  but  whose 
interfacial  area  does  not.  We  shall  refer  to  this  phenomenon  as  the 


^We  us#  the  following  notation:  grad,  div,  and  A  are  the  gradient,  divergence,  and  L^lecian 
operators;  for  F  «  F(t)  and  f  *  ftx.t),  F*  ■  dF/dt  and  f*  *  df/dt;  vol(  )  and  area(-)  denote  the 
volume  and  area  measures. 


formation  of  a  dendrite  with  null  volume. 

Thus  far  we  have  made  no  assumptions  concerning  the  size  of  the 
interfacial  quantities.  Moreover,  even  though  the  hypotheses  underlying 
(1.3)  are  strong,  the  theory  is  exact  in  the  sense  that  the  underlying 
equations  are  fully  compatible  with  the  first  two  laws  of  thermodynamics. 

Next,  we  return  to  the  general  relations  (1.1),  but  in  situations  for 
which  interfacial  energy  and  entropy  are  small.  We  show  that  to  within 
terms  of  higher  order  in  these  quantities. 


lq]m  «  Lv, 
u»-hfc, 


(1.7) 


with  h  as  defined  in  (1.4).  The  relations  (1.7)  are  central  to  the  modem 
work  on  solidification. 1 

We  discuss  a  model  based  on  the  interface  conditions  (1.7)  in 
conjunction  with  assumptions  of  constant  specific  heats  and  constant 
conductivities.  These  assumptions  lead  to  the  equations 


C^-dlvq,  q  =  -k|gradu  in  B,, 


u  *  -  hfc,  [q]-m  =  Lv  on  I,  (1 .8) 

vm  n  -  0  on  dl. 

The  assumption 

C,-C2  (1.9) 

is  common  in  the  literature;  granted  (1.9),  we  are  able  to  establish  the 
following  growth  conditions  for  (1.8): 

(1)  for  an  isolated  boundary, 


(vol(B2)  ♦  CVum)*  =  0, 

(harea(I)  -  umvol(B2)+ (C/2)|(u-um)2)*  i  0; 

B 

(ii)  for  an  isothermal  boundary, 

(h  area(I)  -  uQvol(B2)  +  (C/2)f(u-u0)2l*  i  0. 

B 


(1.10) 


(1.11) 


Here  V-vol(B);  C-C/L;  um  is  the  mean  value  of  u;  in  Case  (ii)  u0  is 
the  constant  boundary-value  of  u . 

A  standard  model  for  solidification  follows  from  (1.8)  when  the  terms 
CjU*  are  neglected: 


Au  =  0,  q  =  -k^gradu 

in  Bjf 

u  ■  -  hfc,  [qlm  =  Lv 

on  I, 

(1.12) 

vmn  =  0 

on  dl. 

Here,  in  place  of  ( 1 . 1 0)  and  ( 1 . 1 1 ),  we  have  the  conditions: 

(i)  for  an  isolated  boundary, 

vol(Bj)*=0,  area(I)*sO;  (1.13) 

(ii)  for  an  isothermal  boundary, 

u0voKB,)*  ♦  harea(I)'iO.  (1.14) 

We  discuss  consequences  of  the  relations  ( 1 . 1 0)  and  ( 1 . 1 1 )  and  of  the 
relations  (1.13)  and  (1.14). 

Because  of  space  limitations,  I  do  not  discuss  the  related  theory  which 
allows  for  variable  concentration.  As  would  be  expected,  the  theory  and 
results  are,  for  the  most  part,  completely  analogous  to  those  presented 
here. 

I  do  not  list  any  but  the  most  basic  hypotheses  of  smoothness,  since 
such  hypotheses  tend  to  obscure  the  main  ideas,  and  since  at  this  stage  of 
the  theory  they  are  academic:  it  is  not  at  all  clear  what  are  appropriate 
function  spaces  in  which  to  analyze  the  underlying  partial-differential 
equations. 


PART  I.  EQUILIBRIUM  THEORY. 


2.  Assumptions.  Terminology. 

2.1.  Constitutive  assumptions. 

We  consider  a  body  B  consisting  of  two  phases  separated  by  an 
interface  I,  and  write  B,  for  the  subregion  of  B  occupied  by  phase  i 

(1*  1,2).  (Cf.  Figure  2).  We  assume  that: 

(Al)  B  is  a  closed,  regular  region  with 
B  bounded  unless  stated  otherwise; 

(A2)  the  B,  are  closed,  regular  regions  with  disjoint  interiors,  and  with 

B-B,UB2,  I-dB,ndB2, 
so  that  B,  and  B2  partition  B. 

We  consider  bodies  whose  behavior  in  equilibrium  is  described  by  five 
fields: 

E(x),  the  bulk  internal  energy  per  unit  volume, 

S(x),  the  bulk  entropy  per  unit  volume; 

T(x),  the  temperature; 

e(x),  the  Interfacial  energy  per  unit  area; 

s(x),  the  interfacial  entropy  per  unit  area. 

We  assume  that  in  each  of  the  phase  regions  Bj  the  material  is  governed 

by  constitutive  equations 

E  =  E|(T),  S  *  S,(T),  (2.1) 

with1 

S,(T)  *  -Fj'(T),  F,(T)  =  E,(T)  -  T  S,(T),  (2.2) 

Fj(T)  being  the  free  energy  of  phase  i.  We  assume  further  that: 

(A3)  the  specif  ic  heats 

C,(T)  -  Ej'(T)  (2.3) 

are  strictly  positive; 

1The  prime  denotes  differentiation  with  respect  to  temperature. 


(A4)  the  free  energies  F,(T)  coincide  at  a  unique  temperature  T  -  TM 
called  the  transition  temperature: 

F/V  •  F2(Tn).  (2.4) 

The  latent  heat  L  is  the  difference  in  energy  between  phases  at  the 
transition  temperature: 

L  -  E2(T„)  -  E,(Tm),  (2.5) 

and  a  third  assumption  we  shall  make  is  that  L  *  0;  by  labelling  phases 
so  that 

phase  2  has  higher  (internal)  energy  at  the  transition  temperature, 

(2.6) 


we  may,  without  loss  in  generality,  write  this  assumption  as: 

(A5)  L  >  0.  (2.7) 

Then,  by  (2.2)2,  (2.4),  and  (2.5), 

S2(  V  ‘  V  V  s  L/ V  °-  (2.8) 

In  view  of  the  agreement  (2.6),  we  shall  adopt  the  following 

Convention.  We  will  refer  to  phase  1  as  the  solid  phase,  phase  2  as  the 

liquid  phase. 

We  do  this  for  convenience  only;  the  theory  is  equally  valid  for 
liquid- liquid  systems,  liquid-vapor  systems,  etc. 

The  next  result,  a  direct  consequence  of  the  our  assumptions,  shows  that 
the  bulk  free-energies  have  the  form  shown  in  Figure  3. 

Properties  of  the  bulk  free-energies. 

(1)  The  bulk  free-energies  are  strictly  concave:  for  T  *  T0, 

F,(T)  <  F,(T0)  -  F,'(T0)(T  -TQ).  (2.9) 

(ii)  The  bulk  free-energies  cross  at  T*Tn;  in  fact, 

F,(T)<F2(T>  for  T  <  T„;  F,(T)>F2(T)  for  T  >  T„.  (2.10) 


-9- 


Figure  3.  Free-energtes  F,(T)  as 
the  transition  temDerat 


mu  Let  u  -  i  -  iM.  men,  as  u  -  0, 


F,(T)-F2(T)*Lu/Tn  ♦  0(u2). 


(2.11) 


Proof.  By  (2.2)  and  (2.3), 

Fj"(T)  -  -  C,(T)/T  <  0, 

which  implies  (i).  Further,  (2.2),  and  (2.8)  imply  that 
F,-(TM)-F2-(TM)-L/Tn>0. 

which  with  (A4)  yield  (ii)  and  (iii).  □ 

The  last  result  shows  that  the  solid  phase  has  lower  free-energy  at 
temperatures  below  transition,  the  liquid  phase  has  lower  free-energy  at 
temperatures  above  transition.  Even  so,  it  is  possible  to  have  the  liquid 
phase  at  temperatures  T  <  TM,  or  the  solid  phase  at  temperatures  T  >  TM; 

in  the  former  case  the  liquid  is  supercooled;  in  the  latter  the  solid  is 
superheated.1  Generally,  we  will  use  the  term  superthermal  to  designate 
either  of  these  two  situations.  Precisely:  phase  i  is  superthermal  at 
temperature  T  if,  for  k  the  other  phase, 

Fk(T)  <  F.(T). 

We  write  constitutive  equations  for  I  analogous  to  those  for  the 
individual  phases.  In  particular,  we  assume  that 


e  -  e(T),  s  -  s(T), 

s(T)  -  -f'(T),  f(T)  -  e(T)  -  Ts(T), 


(2.12) 


(2.13) 


and  that  the  interfacial  free-energy  f(T)  and  the  interfacial 
specific-heat 


c(T)  -  e'(T) 


obey  the  following  hypothesis  for  all  T  >  0: 
(A6)  f(T)  >  0,  c(T)  2  0. 

Then,  arguing  as  above,  we  see  that 

'Cf.  the  interesting  discussion  by  Delves  ( 1974],  Section  3.2. 
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(2. 1 4) 


(2,15) 


m 


f  is  concave 


(2.16) 


(although  not  necessarily  strictly  concave). 

it  is  convenient  to  introduce  the  Gibbs  functions: 

Wj(T,T0)  -  Ej(T)  -  T0Sj(T),  w(T,T0)  -  e(T)  -  T0s(T).  (2. 1 7) 

Clearly, 

F,<T0)  -  Wf(TQ  ,T0).  f  (T0)  «  W(T0,T0).  (2. 1 8) 

In  addition,  we  have  the  following  less  trivial 

Properties  of  the  Gibbs  functions. 

(i )  For  T-  T0, 

F,(T0)  <  wi<T.V*  0  <  f  <V  i  wa.V  (2. 1 9) 

(ii)  For  each  fixed  T0,  each  of  the  Gibbs  functions  (2. 1 7):  has  a  minimum 

at  T  =  T0;  is  monotone  decreasing  on  (0,  T0);  is  monotone  increasing  on 
(T0l«>).  Moreover,  for  the  minimum  and  monotonicity  are  strict 

Proof  We  will  establish  (1)  and  (ii)  for  W,  only;  the  proof  for  w  is 
strictly  analogous.  The  inequality  (2.19),  follows  from  (2.2)  and  (2.9) 

(with  the  roles  of  T  and  T0  reversed).  Also,  by  (2.2),  (2.3),  and  (2.17), 

(d/dDWjdJo)  =  (T  -  T0)Ct(T)/T, 

which  yields  (ii)  for  W..  □ 

Another  consequence  of  the  strict  positivity  of  the  bulk  specific-heats, 
obtained  by  combining  (2.2)  and  (2.9),  is  the  inequality 

E,(T)  -  E,(T0)  >  T0(S,(T)  -  S,(T0)1  (2.20) 

for  T  *  T0;  (2.20)  expresses  the  strict  convexity  of  the  bulk  energies  as 
functions  of  the  corresponding  bulk  entropies. 
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2.2.  Frequently  used  moduli. 

In  this  section  we  list  constitutive  moduli  that  are  used  repeatedly 
throughout  the  paper. 

Modulus 

Section 

F*(T0)  *  bulk  free-energy  of  the  bounded  phase 
minus  the  bulk  free-energy  of  the 
unbounded  phase 

3.2 

b  =  -  f(T0)/F*(T0) 

3.2 

h  •  TMf(TM)/L 

8.1 

a  =  Tms/L 

8.1 

U  *  TMf(T0)/L 

8.1 

13  =  e/L 

8.1 

A  =  h(C,  -  C2)  -  e 

10.1 
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3.  stability. 

3.1.  I  sothermal  boundary. 

Given  a  state  of  B,  that  is,  a  temperature  field  T(x)  and  a  partition 
(B1fB2)  of  B  with  interface  I-dBtndB2,  the  quantities1 


E(T,B},B2)  =  2|Ej(T)  ♦  Jed) , 

B|  1 

S(T,B,,B2)  =  Z/Sjtt)  ♦  Js(T) , 

B;  I 


(3.1) 


represent,  respectively,  the  corresponding  total  energy  and  total 
entropy. 


Remark.  We  will  consider  as  the  underlying  class  of  temperature  fields 
all  functions  T(x)  >  0  on  B  with  the  understanding  that  the  total  energy 
is  and  the  total  entropy  -®°  if  any  of  the  integrals  in  (3.1)  fail  to 
exist.  Further,  many  of  our  statements  concerning  such  fields  will  be 
modulo  a  rearrangement  on  a  set  of  measure  zero;  the  measure,  either 
volume  or  interfacial  area,  will  be  clear  from  the  context. 

We  assume  (for  the  remainder  of  the  section)  that 


the  boundary  is  held  at  (constant)  temperature  TQ. 


The  functional 


F(T,B1,B2)  =  E(T,B,,B2)  -  T0S(T,B},B2)  (3.2) 

then  represents  the  total  free-energy  that  would  be  "recorded"  by  an 
experimenter  not  having  access  to  the  detailed  temperature  distribution  in 
the  interior  of  B.  We  define  as  stable  those  states  of  B  that 


minimize  F(T,B1,B2) .  (3.3) 

Combining  (3. 1 )  and  (3.2)  we  arrive  at  integrals  over  B,  and  I  with 
integrands  equal  to  the  Gibbs  functions  Wj(T,T0)  and  w(T,T0), 

respectively.  Thus,  by  (2. 1 8)  and  (2. 1 9),  a  necessary  and  sufficient 
condition  for  a  minimizer  is  that  T  =  T0,  that  I  =  0,  and  that  B  be 

occupied  by  the  phase  with  lower  free-energy.  Thus,  by  (2. 10)  and  the. 
convention  introduced  in  Section  2,  we  have  the  following,  physically 
obvious, 

1  Here  and  in  what  follows,  X  denotes  summation  over  i  *  1.2. 
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Transition  Theorem.  Assume  that  dB  is  held  at  the  constant 
temperature  T0.  Then  the  stable  states  are  single  phase  with 

temperature  everywhere  constant  and  equal  to  T0.  Moreover:  if  T0  <  T, 
the  solid  phase  is  the  stable  phase;  if  T0>  TM,  the  liquid  phase  is  the 
stable  phase;  if  T0  =  TM,  then  both  phases  are  stable. 


vv 


3.2.  Unbounded  body  under  Isothermal  conditions.  Interracial 
Instabilities. 

As  we  shall  show  In  this  section,  the  theory  Is  more  Interesting  when 
6  Is  unbounded,  In  particular,  we  assume  that: 

B  is  all  of  space. 

We  are,  in  effect,  modelling  a  large  region  whose  boundary  is  held  at  the 
constant  temperature  T0,  and  it  would  therefore  seem  appropriate  to 

define  stable  states  through  (3.3).  The  problem  with  this  definition  Is 
that,  since  B  Is  unbounded,  the  integrals  appearing  In  (3.1)  will  generally 
fail  to  exist.  However,  for  any  partition  (B,,B2)  of  B  and  temperature 

field  T(x)  we  can  consider,  in  place  of  F(T,B1,B2),  the  approximate 

free-energy  FQ(T,B1lB2)  obtained  by  replacing  B,  and  I  by  their 

intersections  with  a  large  ball  0.  Then,  in  view  of  the  discussion  leading 
to  the  Transition  Theorem,  for  any  such  0, 

FQ(To,B1,B2)  <  F0(T,B1,B2) 

for  T  not  identically  equal  to  T0. 

Thus  we  are  led  to  consider  the  problem  (3.3)  with  the  body 
Isothermal  at  temperature  T0.  in  this  case  (3.2)  takes  the  form 

F(T0,B1#B2)  «  2JF,(T0)  ♦  ff(T0).  (3.4) 

We  still  have  the  problem  of  nonexistent  integrals,  but  this  can  be 
circumvented  provided  we  limit  our  discussion  to  situations  in  which 

one  of  the  phase  regions  is  bounded, 

and  this  we  shall  do.  Then,  by  subtracting  the  same  constant  K  from  both 
bulk  free-energies  we  can  give  the  problem  (3.3)  meaning;  indeed,  if  we 
take  K  to  be  the  bulk  free-energy  of  the  unbounded  phase  at  T0  and 

normalize  the  free  energies  in  this  manner,  we  are  led  to  consider,  in 
place  of  (3.4),  the  functional 

F(D)  =  F*(T0)vol(D)  ♦  f(T0)area(dD), 

where  D  is  the  bounded  phase  region,  while  F*(T0)  is  the  bulk 

free-energy  of  the  bounded  phase  minus  the  bulk  free-energy  of  the 
unbounded  phase,  both  evaluated  at  T0.  In  view  of  the  discussion  above, 
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we  define  as  stable  those  bounded  phase  regions  0  that 

minimize  F(D).  (3.5) 

Problem  (3.5)  is  trivial  when  F*(T0)  1 0:  the  solution  Is  D  ■  0.  Thus 
assume: 

F*(T0)  <  0. 

Let  b>0  be  defined  by 

b  -  -  f(T0)/F*(T0),  (3.6) 

so  that  (3.5)  has  the  form 

minimize  -vol(D)  ♦  b  area(dD):  (3.7) 

By  considering  a  sequence  of  balls  with  radii  tending  to  infinity,  it  is 
clear  that  the  infimum  of  the  functional  in  (3.7)  is  and  there  are  no 
stable  phase  regions.  However,  even  though  (3.7)  has  no  solution,  we  can 
study  the  types  of  Instabilities  compatible  with  our  model  by  studying 
minimizing  sequences;  that  is,  sequences  (Dn)  on  which  the  functional 

in  (3.7)  tends  to  -«>.  Let 

v„  ■  volCD')'  a„  *  areaCdDJ. 

n  n  n  n 

Then  (Dn)  is  a  minimizing  sequence  if  and  only  if  (Dn)  is  consistent  with 

the  area-volume  limit: 

v„  -ba„-*r<»  as  n  **■«>. 

n  n 

Since  b  >  0,  an  immediate  consequence  of  this  condition  is 
vn-*«>  as  n -*•«>. 
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Remark.  By  (2. 1 1 ),  the  modulus  b  has  the  following  approximate  form 
for  u0  -  T0  -  Tm  small: 

b  *  T„f(T0)/Llu0l.  (3.8) 

Mullins  and  Sekerka  [1963]  estimate  that  f(Tn )/L  has  10"7cm.  as 
order-of -magnitude.  Taking  this  as  an  estimate  of  f(T0)/L,  and  using  the 
value  10'2  given  by  Delves  ([1974],  p.  80)  as  a  "small  value"  of  IuqI/T^, 

we  are  led  to  the  estimate  of  10"5  cm.  as  a  reasonable  "large  value"  of  b. 
(Of  course,  b  can  be  as  large  as  we  wish;  we  simply  take  T0  sufficiently 

close  to  1^.) 

The  area-volume  limit  restricts  the  growth-rate  of  interfacial  area, 
but  even  so  it  is  possible  to  construct  minimizing  sequences  in  which  the 
interface  exhibits  interesting  behavior.  We  now  give  some  examples  (cf. 
Figure  4). 

1.  Dendrite  spike. 

Let  Dn  be  a  right  circular  cylinder  of  radius  r  and  height  n,  with 
r  independent  of  n  and  r  >  2b  (cf.  the  Remark).  Then 

vn  -  ban  =  0(n)  (3.9) 


and  the  area-volume  limit  is  satisfied.  Thus  (Dn)  is  a  minimizing 
sequence.  We  could  let  r  ■  rn  depend  on  n,  as  long  as  inf  rn  >  2b. 


2.  Dc.iurite  star. 

Let  Dn  be  the  union  of  N(n)  spikes,  all  of  height  n  and  radius  r>2b, 


where  N(n)  is  finite,  but  arbitrary.  Then  (3.9)  remains  valid  and  (Dn)  is  a 


minimizing  sequence. 


3.  Simple  dendrite  tree. 

Let  Dn  consist  of  N(n)  spikes,  each  of  radius  r,  but  with  varying 
heights  which  total  l(n),  where  l(n)-*-oo  as  r Then,  for  large  n, 


vn  *  TTf^Kn)  +  0(1),  an  =  2TTrl(n)  +  0(1), 


and  the  area-volume  limit  is  again  satisfied. 
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v  v.  /,  /, 


ii 


simple  dendrite  tree 


,1 

$ 

'v? 


o(n2)  spikes  of  length  n 
covering  a  sphere  of 
radius  n 


prickly  ball 


Figure  4.  Some  examples  of  minimizing  sequences  for  (3.5)  when  the 
unbounded  phase  is  superthermal. 
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4.  Dendrite  tree  with  branches  of  varying  radii. 

Let  Dn  consist  of  a  stalk,  a  dendrite  spike  of  radius  rft  and  height  n, 

and,  growing  from  it,  np  spikes  arranged  as  follows:  for  each  m, 
m  *  1 ,2,...,n,  there  are 

p  spikes  of  radius  rn/m  and  height  n/m. 

Then 

vn  *  ,hrn2hp(p"1  +  2~3  ♦  3"3  ♦  ...  ♦  n-3)  +  0(1), 
an  =  2Tirnnp(p_1  *  2'2  ♦  3'2  +  ...  ♦  n-2)  +  0(1 ). 

For  k  -  2,3, 

OB 

*  p"1  +  2  m-* 
m*2 

is  well  defined.  Let  rn-*r,  with  r  finite  or  +°°.  Then,  if 

rAj/Aj  >  2b, 

the  area-volume  limit  is  satisfied. 

5.  Prickly  ball. 

Here  Dn  consists  of  a  sphere  of  radius  n  covered  with  o(n2)  spikes  of 
length  n  and  radius  r.  Then, 

v-  ■  (4/3)Ttn3  +  o(n3),  an  *  o(n3) 

and  the  area-volume  limit  is  satisfied. 

The  sequences  described  above  furnish  simple  examples  of  minimizing 
sequences  and  display  the  richness  of  the  underlying  theory;  less  trivial 
examples  which  more  closely  resemble  dendritic  behavior  are  easily 
constructed. 

To  state  the  results  of  this  section  succintly,  recall  that,  in  the 
terminology  of  Section  2, 

F*(T0)  <  0  if  and  only  if  the  unbounded  phase  is  superthermal; 

F*(T0)  >0  if  and  only  if  the  bounded  phase  is  superthermal. 


Theorem.  Assume  that  the  body  is  isothermal  at  temperature  T0 . 

( i)  Consider  (3.5)  with  the  unbounded  phase  superthermal  Then  (3.5)  has 
no  solution)  A  sequence  (Dn)  is  a  minimizing  sequence  if  and  only  if  it 

obeys  the  area-volume  limit  and  therefore  has  corresponding  volumes 
tending  to  infinity  These  conditions  allow  for  interfacial  instabilities 
such  as  those  exhibited  by  the  sequences  described  above.. 

( i  i )  Consider  (3.5)  with  T0  =  TH,  or  with  the  bounded  phase  superthermal. 

Then  the  only  solution  is  D  =  0. 


Remarks. 

(1)  It  would  be  interesting  if  one  could,  in  some  sense,  geometrically 
classify  the  sequences  consistent  with  the  area-volume  limit. 

(ii)  By  definition,  a  phase  is  superthermal  if  it  has  higher  free-energy  at 
the  underlying  temperature.  In  this  sense  the  superthermal  phase  is 
unstable,  and,  as  the  theorem  shows,-  the  other  phase  ultimately  occupies 
the  entire  body. 

(iii)  Mullins  and  Sekerka  [1963],  working  within  the  dynamical  theory 
discussed  in  Section  1 1,  have  established  the  instability  of  the  interface 
for  (infinitesimal  perturbations)  of  a  sphere  solidifying  in  a  supercooled 
melt.  Assertion  (i)  of  the  theorem,  with  [Dnl  the  prickly  ball,  gives  an 

analog  of  the  Mullins-Sekerka  instability  within  the  framework  of  the 
equilibrium  theory. 

(iv)  The  Transition  Theorem  and  the  theorem  above  show  exactly  what  is 
needed  for  the  interface  to  be  unstable,  namely: 

(a)  an  infinite  region,  for  otherwise  the  boundary  captures  the  instability; 

(b)  the  unbounded  region  in  a  superthermal  state,  so  that  it  is  enveloped  by 
the  more-stable  bounded  region. 


^  Actually,  in  this  case  the  empty  set  is  a  local  minimum  of  (3.5),  since  (3.5)  is  positive  for  all  D  of 
sufficiently  small  volume. 
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3.3.  Isolated  body. 

As  a  consequence  of  the  transition  theorem,  stable  states  are 
necessarily  single  phase  when  the  boundary  is  held  at  constant 
temperature.  We  now  show  that  coexistent  phases  are  possible  when, 
instead,  the  body  is  isolated.  For  convenience,  we  return  to  our  original 
assumption  that  B  be  bounded. 

For  an  isolated  body  the  total  energy  is  fixed.  Writing  E0  for  the 

value  of  this  energy,  the  corresponding  stable  states  are,  by  definition, 
those  that 

maximize  S(T,B,,B2)  subject  to  E(T,B,,B2)  =  EQ.  .  (3.10) 

To  solve  the  problem  (3.10)  we  assume  that  the  bulk  and  Interfacial 
specific  heats  are  strictly  positive,  and  that 

Ej(0+)  is  finite,  Ejf00) =  +°°;  ^(0+)  =  -<*.  (3.1 1) 

Also,  to  insure  that  the  constraint  equation  in  (3.10)  has  a  solution,  we 
assume  that 

EQ  >  EjtO-OV  for  at  least  one  value  of  i.  (3.12) 

Our  next  result  shows  that  (3. 1 0)  has  a  solution;  in  the  statement  and 
proof  of  this  theorem  we  use  terminology  introduced  in  the  Appendix. 
Also,  for  convenience,  we  write 

JM=  (EjCT^V,  E2(Tm)V).  (3.13) 

Theorem.  Consider  an  isolated  body.  Then  there  exist  stable  states 
corresponding  to  any  E0  consistent  with  (3.12).  Moreover: 

(1)  Each  stable  state  has  constant  temperature. 

(ID  Each  two-phase  stable  state  has  minima /  interface;  hence  the 
interface  is  a  surface  of  constant  mean  curvature  which  meets  dB 
orthogonally 

(ill)  For  f  (TM )  sufficiently  small  there  is  a  nonempty  interval 
JcJM  such  that  for  E^j  all  stable  states  are  two  phase. 

Proof.  Our  proof  will  proceed  in  three  steps. 

Step  t:  Proof  of  (i).  Fix  B,,  B2  (and  hence  I),  let 

v^voKBj),  a  =  area(I), 
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and  consider  (3.10)  as  a  problem  for  T.  Let  us  agree  to  call  a  temperature 
field  T  admissible  If  T  is  consistent  with  the  energy  constraint  In 
(3.10);  phasewise  constant  If  T  is  constant  on  B,,  constant  on  B2,  and 

constant  on  I.  Further,  for  each  temperature  field  T  define  numbers 
Tj  >  0  and  T,  >  0  by 

EttjJv,  =  |E,(T),  e(T,)a  -  Jed). 

Bf  1 

(Because  of  the  strict  monotonicity  of  E,(T)  and  e(T),  this  definition  is 
valid.)  Then,  in  view  of  (2.20), 

f(S,(T)  -  S,(T,))  -  /(S,(T)  -  S((T()  -  T,~*(Ej(T)  -  E^T,)))  <  0 

Bi  Bl 

if  T  Is  not  identically  equal  to  T,  on  B,;  a  similar  argument  applies  to 

s(T)  on  I.  Assume  that  T  is  admissible,  but  not  phasewise  constant. 
Then,  by  the  above  inequalities,  the  phasewise-constant  field  with  values 
T,  on  B,  and  T,  on  I  is  admissible  and  has  total  entropy  strictly  larger 

than  that  for  T.  Thus 


the  set  of  stable  states  Is  contained  in  the  set  of 
states  with  phase  wise-constant  temperature  fields. 


(3.14) 


Continue  to  hold  Bv  Bj  (and  hence  I)  fixed.  Consider  the  problem, 
motivated  by  (3.14),  of  finding  numbers  Tj>0  and  T,>0  that 


maximize  S^TjJv,  ♦  S2(T2)v2+  s(T,)a 
subject  to  E,(T1  )v,  ♦  E2(T2)v2  ♦  e(T,)a  *  E0. 


(3.15) 


In  view  of  (3.1 1)  and  the  continuity  of  the  functions  appearing  in  (3.15), 
this  problem  has  a  solution  and,  using  a  Lagrange-multiplier,  it  is  a  simple 
matter  to  verify  that  solutions  have  T^Tj  *  T,.  Thus 


the  set  of  stable  states  Is  contained  in  the  set 
of  states  with  constant  temperature  fields, 

which  establishes  assertion  (i)  of  the  theorem. 


(3.16) 


Step  2:  Existence  of  stable  states  and  proof  of  (11).  Consider  the  problem 
of  finding  numbers 


s 


$ 

s 


I& 

53 


T  >  0,  0  i  v  i  V,  a  *  amin(v) 


(3.17) 


that  solve  the  problem 

minimize  W(a,v  J  ;T0)  =  W,(T,T0)v  +  W2(T  ,T0)(  1  -v)  ♦  w(T  ,T0)a 


(3.18) 


subject  to  E,(T)v  ♦  E2(T)(  1  -v)  ♦  e(T)a  =  E0, 


where  amin(v)  is  the  minimal-area  function  defined  in  the  Appendix,  W, 
and  w  are  the  Gibbs  functions  (2.17),  and  T0  >  0  is  arbitrary,  but  fixed. 
Because  of  (2.17)  and  (3.16),  every  stable  state  (T.B^Bj)  yields  a  solution 
(a,v,T)  of  (3.18)  with 


vsvol(Bt),  a  =  area(l); 


(3.19) 


conversely,  every  solution  (a,v,T)  of  (3. 1 8)  generates  -  in  the  obvious 
sense  -  a  stable  state  provided  we  can  find  a  partition  of  B  consistent 
with  (3.19).  But  by  Theorem  1  of  the  Appendix,  such  a  partition  can  always 
be  found  if 


3  “  amin(v>- 


(3.20) 


Thus  to  complete  Step  2,  It  suffices  to  show  that  (3.18)  has  a  solution,  and 
that  every  solution  obeys  (3.20). 

Because  of  (2. 1 7)  and  (3. 1 1 ),  w(T,T0)  is  bounded  below,  while 

W,(0*,T0)  =  +~,  W,(+«,T0)  =  ♦<*>; 

thus  (3.18)  has  a  solution.  Let  (a,v,T)  be  such  a  solution.  Then  (a,v,T)  is 
a  solution  for  any  T0  >  0.  Choose  T0  >  T.  Assume  that  a  >  amjn(v).  Then, 

because  of  the  strict  monotonicity  of  the  Internal  energies,  we  can  find 
numbers  a*  and  T*  such  that  T  <  T*  <  T0  and  amjn(v)  <  a*  <  a;  and  such 

that  (a*,v,T*)  satisfies  the  constraint  equation  in  (3.18).  Further, 
because  of  (ii)  of  the  theorem  containing  (2.19), 

W(a*,v,T*;T0)  <  W(a,v.T;T0), 

which  contradicts  our  assumption  that  (a,v,T)  solve  (3. 18).  Thus  (a,v,T) 
obeys  (3.20). 


nl 

1 
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Step!  Proof  of  (iii).  Let 

WM  -  W/W  -  W2(  Tm.T„),  am„  -  sup  (a„Jv>:  0 !V£V), 
and  let  T,  »  T,  (E0)  be  the  unique  solution  of 


E,(T,)V  -  E0: 


Then,  as  is  clear  from  Figure  5  and  the  properties  of  the  Gibbs  functions, 
there  is  a  nonempty  interval  J<r  J„  such  that  for  EQeJ, 


W<(T|(W>  W^f(TM)amax  >  W(amin(v),v,TM;TM). 


(3.21) 


The  left  side  of  this  inequality  is  the  total  Gibbs-energy  (3.18),  (with  T0  ■ 
Tm)  for  the  two  admissible  single-phase  states.  On  the  other  hand,  for 

0  <  v  <  V  the  right  side  is  the  total  Gibbs-energy  for  a  (not  necessarily 
admissible)  two-phase  state  at  temperature  T  *  TM.  Thus,  if  we  can  find 

a  v€(0,V)  that  satisfies  the  energy  constraint  (3.1 8)2  with  T  *  TM,  we 

will  have  exhibited  an  admissible  two-phase  state  with  lower  total 
Gibbs-energy  than  the  two  admissible  single-phase  states,  so  that  the 
solution  of  (3. 1 8)  is  necessarily  two  phase.  Thus  to  complete  the  proof  we 
have  only  to  find  a  ve(0,V)  such  that 


E,(TM)vtE2(TM)O-v)  +  e(T)amin(v)  =  E0. 


(3.22) 


But  since  E0£jM,  and  since  the  left  side  of  (3.22)  varies  continuously 
from  E2(Tm)V  at  v  a  0  to  E,(TM)V  at  v  =  V,  we  may  use  (3.13)  to 
conclude  that  such  a  v  exists.  □ 


|tv«] 


PART  II.  GENERAL  THEORY. 


4.  Preliminaries. 

We  now  allow  the  interface  I(t)  and  (hence  also)  the  phase  regions 
B,(t)  to  depend  on  the  time  t.  We  use  the  following  notation  (Figure  2): 

n(x)  is  the  outward  unit  normal  to  dB;  m(x,t)  is  the  outward  unit  normal 
to  Bt(t)  on  I(t);  v(x,t)  is  the  normal  velocity  of  the  Interface  in  the 

direction  m(x,t);  Mx,t)  Is  the  sum  of  principal  curvatures  of  I(t) 
with  fc>0  when  the  center  of  curvature  lies  toward  B,;  cot(m.n)  is  the 
cotangent  of  the  angle  between  m  and  n: 

cot(mfn)-p/(1-p2)1/2,  p-mn.  (4.1) 

For  any  field  g(x,t),  let  gc0  denote  the  limit  of  g  as  the  Interface  Is 
approached  from  points  of  B,,  and  let 

tgl  -  9(2)  -  9(i)  (4.2) 

denote  the  jump  in  g  across  I.  We  then  have  the  standard  identity 

( Jg  )•  •  fg*  ♦  Jg*  -  JCglv.  (4.3) 

6  a,  &2  1 

A  less  trivial  result  holds  for  integrals  over  the  interface.  Let  g(x,t) 
be  defined  for  each  x  in  I(t)  and  all  t.  Choose  a  time  B  and  a  point  y 
in  1(G),  and  let  x(t)  denote  the  curve  that  passes  through  y  at  time  J3 
and  has 

x*(t)  -  v(x(t),t)m(x(t),t) 
for  all  t.  Then 

gA(y,J3)  *  (d/dt)g(x(t),t)|l-B  (4.4) 

represents  the  time-derivative  of  g  following  I. 


Transport  Theorem.1 


( jg  >*  =  jgfcv  ♦  |gA  -  Jjgvcot(mn). 


(4.5) 


The  first  term  on  the  right  side  of  (4.5)  represents  a  change  in  the 
integral  of  g  due  to  the  curvature  of  the  interface,  the  second  term  gives 
the  change  due  to  changes  in  g,  the  third  term  represents  a  flux  of  g 
across  the  boundary  of  B. 

An  immediate  consequence  of  (4.5)  is  the  transport  theorem  for  area: 


area(I)*  *  ffcv  -  Jvcot(m-n). 
I  at 


(4.6) 


Remark.  By  (4.1),  cot(m‘n)  =  «  when  msuv  Thus  for  (4.5)  and  (4.6)  to 
be  valid  it  is  necessary  that 


on  ai:  v  *  0  when  m  *  in. 


(4.7) 


and  that  this  limiting  value  of  v  be  approached  sufficiently  rapidly  from 
the  interior  of  I.  This  condition  will  be  a  tacit  asumptlon  in  what 
follows. 


Let  S  denote  an  arbitrary  subsurface  of  dB,  let  r  >  0,  and  let 


Dr  =  (x:  x  =  y  -  J3n(y),  ye  S,  0  i  B  i  r). 


Sr  -  (x:  x  -  y  -  rn(y),  yeS). 


Then  Dr  and  Sr,  respectively,  are  the  region  and  surface  obtained  by 
r-transporting  S  along  the  normal.  Here  we  will  always  write  n,. 
for  the  outward  unit  normal  to  Dr  on  Sr  (Figure  6). 


'Williams  ( 19051.  In  this  connection  see  also  Scrtven  [I960),  noeckel  [  1975). 
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Figure  6.  The  region  Dr  and  surface  Sp  obtained  by  r- transporting  s 
along  the  normal. 


5.  First  two  laws.  Additional  constitutive  assumptions. 

To  discuss  the  thermodynamics  of  B  we  allow  the  functions  E(x,t), 
S(x,t),  T(x,t),  e(x,t),  and  s(x,t)  to  depend  on  the  time  t,  and  we  introduce  a 
new  field:  q(x,t),  the  heat  flux  per  unit  area. 

Let  D  be  an  arbitrary  fixed  subregion  of  B,  and  let 

ID  *  In  D;  Cp«dI0\aB,  (5.1) 

so  that  dID  is  the  closed  curve  which  marks  the  intersection  of  I  with 
dD,  while  Cp  is  the  portion  of  this  curve  which  lies  in  the  interior  of  B. 
The  internal  energy  of  D  is  given  by 

JE  +  Je , 

D  «d 

Also,  because  of  the  (possible)  motion  of  the  interface  relative  to  D, 
there  is  a  flux  of  interfacial  energy  across  the  boundary  of  D  (and  out  of 
D)  of  amount 

/evcot(m  n) , 

CD 

where  n  is  the  outward  unit  normal  to  dD,  while  cot(m  n)  is  defined  in 
(41)  (cf.  Figure  7).  We  integrate  over  Cp,  rather  than  ai0,  since  there  is 

no  loss  of  interfacial  energy  across  the  boundary  of  B .  Finally,  there  is  a 
flux  of  heat  into  D  of  amount 

-  Jq  n. 

dD 

In  view  of  the  above  discussion,  the  first  law  for  D  takes  the  form: 

( fE  +  Je  )•  +  Jevcot(m  n)  *  -  Jq  n.  (5.2)1,2 

D  Ip  Cp  dD 

Similarly,  D  has  internal  entropy  and  there  is  a  flux  of  interfacial 
entropy  across  the  boundary  of  D.  Also,  the  flow  of  heat  is  accompanied 
by  a  flux  of  entropy  of  amount 


1  Fernandez -Diaz  and  Williams  [  1979).  An  earlier  version  is  contained  in  the  work  of  and  Moeckel 
( 19751,  but  Moeckel  introduces  a  somewhat  artificial  material  description  of  the  interface.  Other 
discussions  of  the  first  law  are  given  by  Fisher  and  leitman  1 1968),  Murdoch  [  1976),  Wollkind  (1979), 
and  Rogers  (1983). 

^Surface  tension,  in  our  theory  equal  to  the  interfacial  free-energy.  does  not  enter  the  first  law; 
indeed,  since  our  model  is  rigid,  bulk  and  interface  forces  do  no  work. 

-30- 


I 


We  therefore  write  the  second  law  for  D  as 


( fs  ♦  fs)*  +  fsvcot(mn)  >  -jT_1qn. 

6  fD  CD  30 


(5.3) 


The  term 


({5  ♦  fs)*  ♦  Jsvcot(mn)  ♦  jT_,qn 

6  f0  c0  30 


(5.4) 


represents  the  entropy  production  in  D,  and  the  second  law  (5.3)  Is 
simply  the  requirement  that  P(D)  i  0. 

In  addition  to  the  constitutive  assumptions  introduced  in  Section  2.1, 
we  assume  that  the  heat  flux  q  is  governed  by  the  constitutive  equation 


q  =  -kj(T)gradT 


(5.5) 


in  each  phase  region  B<,  with 


(A7)  conductivities  k4(T)  >  0. 


Concerning  the  interface  I,  we  assume  that 
(A8)  T  is  continuous  across  the  interface, 


(5.6) 


so  that  the  interfacial  temperature  is  defined  unambiguously.  In  addition, 
we  assume  that  the  interface  produce  no  entropy: 

(A9)  if  the  volume  of  0  tends  to  zero,  then  P(D)  tends  to  zero.  (5.7) 


6.  Partial-differential  equations.  Jump  conditions. 

In  this  section  we  derive  consequences  of  the  requirement  that  the 
first  two  laws  (5.2)  and  (5.3)  be  satisfied  in  every  subregion  D  of  S 
We  assume  that  the  fields  E,  S,  and  T  are  smooth  away  from  the 
Interface  I.  Then  (5.2)  applied  to  an  arbitrary  region  D  lying  in  B,,  in 

conjunction  with  the  constitutive  equations  (2.1)  and  (5.5),  lead  to  the 
classical  partial-differential  equations 

E}(T)* -  -divq,  q  -  -kj(T)gradT,  (6. 1 ) 

or  equivalently,  by  (2.3), 

C,(T)  T* »  div(k,(T)gradT)  (6.2) 

In  B,  for  all  time. 

Our  next  step  is  to  derive  appropriate  jump  conditions  across  the 
interface.  Although  the  temperature  is  continuous,  the  fields  E,  S,  and  q 
will  generally  suffer  jump  discontinuities  across  the  interface.  We  also 
•allow  q(x,t)  to  become  infinite  as  x  approaches  dl(t);  in  fact,  we 
assume  the  existence  of  a  scalar  function  Q(x,t)  on  dl(t)  such  that: 
given  any  subsurface  S  of  dB,  if  Sp  is  obtained  from  S  by 

r-transporting  S  along  the  normal,  then 


Jqnp  Jqn  ♦  JQ ,  jT-,q-np  jr’q* n  ♦  JT-,Q 

sp  s  snai  s,,  s  snai 

as  r-»0. 

The  next  result  is  central  to  our  investigation. 
Interface  Theorem.  The  Interfaclal  balance  laws 
-v[E]  ♦  vfce  *  eA  *  -[q]m, 

-vis]  ♦  vfcs  ♦  sA  ■  -r’lqj-m 
hold  on  I  for  all  time.  The  contact  condition 
vm-n  ■  0. 

holds  on  ai  Tor  an  time. 


(6.3) 


(6.4)1 


(6.5) 


(6.6) 


Proof.  Note  first  that,  by  (6. 1 )  and  the  identity  E.(T)#  =  TS^T)*, 


'Cf.  Moeckel  ( 1975],  Fernandez -Oiaz  and  Williams  1 1979). 


-jj- 


5j(T)*  =  -  div(T-1q)  ♦  r,  r  =  -T'2q*gradT 


(6.7) 


in  Bj  for  all  time. 


Next,  the  laws  (5.2)  and  (5.4),  for  an  arbitrary  subregion  D  of  B, 
have  the  form 


( }G  ♦  jg  Y  +  Jgvcot(mn)  =  -jpn  ♦  R(D), 


(6.8) 


with  R-0  for  (5.2)  and  R-P  for  (5.4),  so  that  in  either  case  R  is 
consistent  with  (5.7).  By  (4.3)  and  (4.5),  the  left  side  of  (6.8)  is  equal  to 


JG*  ♦  JG*  ♦  Jtgfev  +  gA  -  [G]v}  -  Jgvcot(mn), 

6 ,  Do  In  Sn 


(6.9) 


where  Dj  is  the  portion  of  D  in  Bj,  while 


Sp  =  dlpfi  dB 


is  the  portion  of  dip  in  contact  with  dB. 

Assume  first  that  dD  does  not  touch  dB.  Then  we  may  use  (6.1 ),  (6.7), 
and  the  divergence  theorem  to  conclude  that 


jG*  =  -jpn  +  Z(Dj), 


(6.10) 


where  Z(DS)  tends  to  zero  with  the  volume  of  Dj.  Further, 


J[p]-m  -  ( J  -  |  -  |  )(p  n) , 


(6.11) 


Iq  dD  dD  |  3D2 

and,  since  Sp  is  empty,  if  we  combine  (6.7M6.1 1 )  we  find  that 


j{£,v  ♦  gA  -  [Gjv  ♦  [p] m)  =  R(D)  -  Z(D.)  -  Z(D2). 


(6.12) 


By  (5.7),  if  we  shrink  D  to  the  interface,  the  right  side  of  (6.12)  tends  to 
zero;  since  ID  is  essentially  arbitrary,  this  yields  (6.4)  and  (6.5). 

Now  choose  a. subsurface  S  of  dB,  apply  (6.8)  and  (6.9)  to  the  region 
Dr  obtained  from  S  by  r- transporting  S  along  the  normal,  and  let  r 

approach  zero;  by  (5.7),  (6.3),  and  the  fact  that  S  is  arbitrary,  the  result 


evcot(m-n)  -  Q,  svcot(m-n)  -  T_1Q, 


(6.13) 


so  that,  by  (2,13), 


f(T)vcot(m-n)  =  0, 


"  w  "  rf  •v'y  *  a  0  , M  *  u  M  _  W  <■’■*■**  *  1 " >  ^  *  *  *  *  *  «  *  »  W  *  •  %  0^1  \  *„  %  _■,»**.  f  JF 


which  yields  (6.6),  since  f  >  0.  □ 


Remarks. 

(i)  It  is  clear  from  (6.6)  and  (6.13)  that  Q  -  0,  so  that  (apart  from  a 
possible  jump  discontinuity)  the  heat  flux  is  nonsingular  along  dl. 

(il)  We  could  also  allow  for  surface  energies  ej  and  surface  entropies  s, 

between  the  phases  and  dB.  Were  we  to  do  this,  then,  assuming  e,  and  s, 
are  constant,  (6.6)  would  be  replaced  by  the  more  general  condition 

v(f  ,(T)  -  f2(I)  -  (m-n)f(T)}  *  0,  fs(T)  *  e.  -  Ts, 

on  dl.  Moreover,  in  this  case  the  heat  flux  is  generally  singular  on  dl 
with 

Q  «  -v(eJ  -  e2  -  (m-n)e(T)Kl-  (m-n)2)-1/2. 

As  an  immediate  corollary  of  (6.4)  and  (6.5)  we  have,  for  v  *  0, 
Balance  of  free  energy.  On  1,  for  all  time, 

[F]  8  fcf,  (6.14) 

or  equivalently 

T  *  ([El  -  fce)/([S]  -  fcs).  (6.15) 

Proof.  By  (2.13)  and  (4.4),  eA  s  TsA;  thus  if  we  eliminate  the  heat  flux 
from  (6.4)  and  (6.5),  we  arrive  at  (6.14),  and  solving  this  for  T  yields 
(6.15).  □ 

Remarks. 

(i)  If  e*s*0,  which  is  the  assumption  underlying  the  classical  theory  of 
melting,  or  If  MO,  then,  as  a  direct  consequence  of  (6.15), 

T  *  (E)/[S)  (6.16) 


and 


[F] =  0.  (6.17) 

Thus,  granted  the  assumption  that  the  free  energies  F,(T)  coincide  at 
exactly  one  temperature,  we  find,  as  a  consequence  of  the  laws  of 
thermodynamics,  that  in  the  absence  of  interfacial  structure  the 


free-boundary  condition  (6.15)  must  necessarily  reduce  to  the  classical 
condition 

T  -  Tm  on  I. 

(ii)  Even  though  the  notion  of  force  does  not  enter  our  theory,  we  can 
define  pressure  by  noting  that,  since  F.(T)  is  the  free-energy  per  unit 

volume,  VFj(T)  is  the  free-energy  a  region  of  volume  V  would  have  if 

occupied  solely  by  phase  i;  differentiating  this  expression  with  respect  to 
volume  yields  the  negative  of  the  pressure: 

P,  -  -F,(T). 

A  similar  argument  identifies  f(T)  with  surface  tension  <f.  With  these 
definitions  (6.14)  is  the  Gibbs-Thompson  relation1 

P,  -p^fctf. 

It  must  be  emphasized  that  our  derivation  of  (6.14)  makes  no  use 
whatsoever  of  the  notion  of  equilibrium. 

(iii)  It  is  not  difficult  to  show  that 

P(B)  =  X  kj  jr2lgradTI2, 

Bi 

so  that,  as  would  be  expected,  the  entire  production  of  entropy  is  due  to 
heat  conduction. 

We  close  this  section  by  listing  what  might  be  considered  a  complete 
system  of  field  equations  and  free-boundary  conditions  for  the  general 
nonlinear  theory:  namely,  (6.2),  (6.4),  (6.6),  and  (6.15).  We  do  not  need 
(6.5),  since,  granted  (6.4),  the  conditions  (6.15)  and  (6.5)  are  equivalent. 


The  system  is,  then, 

C,(T)  T*  =  -divq, 

q  =  -k.(T)gradT 

in  B., 

(6.18) 

T  «  (IE)  -  fce)/([S]  -  fcs), 

[q]m  -  vlE]  -  vfee  -  eA 

on  I, 

(6.19) 

vm  n  -  0 

on  dl, 

(6.20) 

with  E,  5,  e,  and  s  given  by  the  constitutive  equations  (2.1),  (2.2),  (2.12), 
and  (2.13). 

’Gibbs  (1878).  Eq.  (500). 


7.  Boundary  conditions.  A  bound  on  Interfaclal  area. 

Throughout  this  paper  we  will  be  concerned  with  the  following 
boundary  conditions: 1 

(I)  isolated  boundary: 
q*n*0  on  dB; 

(II)  isothermal  boundary: 

T  *  T0  on  dB  (T0  constant). 


(7.1) 

(7.2) 


These  boundary  conditions,  when  combined  with  (5.2),  (5.4),  and  (6.6), 
yield  the 

Global  growth  conditions.  if  the  boundary  is  isolated. 


(7.3) 


P(B)  i  0. 


if  the  boundary  is  isothermal  or  isolated, 

{ f(E-T0S)  ♦  J(e-T0s)  Y  i  ~T0P(B)  S  0,  (7.4) 

B  I 

with  the  constant  T0  i  0  arbitrary  when  dB  is  isolated 

Remark.  Note  that  the  variational  characterizations  (3.3)  and  (3.10)  can 
be  deduced  from  the  above  theorem  under  the  assumption  that  each  process 
consistent  with  boundary  conditions  (7. 1 )  or  (7.2)  approach  a  stable 
equilibrium  as  t -*■«>. 

If  we  integrate  (7.4)  from  the  initial  time  t*0  to  an  arbitrary  time, 
we  conclude,  with  the  aid  of  (2. 19),  that 

FjCTqJvoKB,)  ♦  F2(T0)vol(B2)  ♦  f(T0)area(I(t))  i  ft 0), 
ftO)  -  the  initial-value  of  J(E-T0S)  ♦  J(e-T0s). 

B  i 


We  do  not  mean  to  Imply  that  these  are  the  only  boundary  conditions  or  interest.  In  fact,  here  and  In 
subsequent  sections,  all  results  for  an  isothermal  boundary  are  valid  without  change  for  the  boundary 
condition:  T*T0  on  Aj.  q  n  ■  0  on  A2.  with  Af  complementary  (nontrivial) subsets  of  dB. 


Thus,  letting 


Fmin(To) "  V  -  vol(B), 

we  are  led  to  the  following 

Bound  on  interfacial  area.  Assume  that  the  boundary  is  isolated or 
isothermal.  Then  for  all  time  t, 

area(I(t))  <  f(T0)-i(FK»  -  VFm|n(T0)},  (7.5) 

with  T0  arbitrary  when  dB  is  isolated 
Note  that,  in  addition  to  (8.5),  we  have  the  estimate 
VFmin<V  *  J<E-T0S)  i  fXO). 

B 

Intuitively,  one  type  of  dendrite  might  be  a  “set"  with  vanishingly  small 
volume,  but  nonzero  surface  area.  The  following  definition1  makes  this 
idea  precise. 

We  say  that  phase  1  approaches  a  dendrite  with  null  volume  if: 
lim  vol(B,(t))  =  0,  llm  inf  area(I(t))  >  0.  (7.6) 

t-**»  t+«* 


Remark.  The  inequality  (7.5)  asserts  that  area(I(t))  is  bounded 
uniformly  in  t  by  a  constant  which  depends  only  on:  (i)  the  bulk  and 
interfacial  free-energies  at  T0;  (ii)  the  initial  data.  At  first  glance  it 

might  seem  that  this  bound  on  interfacial  area  precludes  the  formation  of 
an  interface  which  is  too  wild.  But  this  is  not  so.  Indeed,  dendrites  with 
null  volume  are  completely  consistent  with  (7.5).  Moreover,  as  we  shall 
see,  the  theory  discussed  here  will  predict  such  instabilities,  but  only  on 
a  very  small  length  scale. 


'a  more  general  definition  would  use  the  limit  t-*t*  with  t#  finite  or  infinite. 


8.  Quasi-static  theory. 

8.1.  Basic  equations.  Growth  conditions. 

In  this  section  we  attempt  to  model  situations  In  which  the  interface 
moves  slowly  compared  with  the  time  scale  for  heat  conduction.1  The 
constitutive  hypothesis  underlying  our  model  is  that  -  in  each  of  the 
phases  -  both  the  bulk  energy  and  the  bulk  entropy  are  constant.  Since  the 
two  phases  coexist  at  the  interface,  the  constant  values  of  energy  and 
entropy  should  be  the  values  in  the  real  material  at  the  transition 
temperature.  Thus,  by  (2.5)  and  (2.8),  and  since  there  is  no  loss  of 
generality  In  allowing  the  energy  and  entropy  to  vanish  in  phase  1 ,  we 
base  our  model  on  the  constitutive  assumptions: 


S,  -  O.  S2  =  L/Tn, 

with  L  >  0  the  latent  heat  (2.5).  In  addition  to  (8. 1 ),  we  suppose  that2 
e,  s,  and  k,  are  strictly-positive  constants;  f(TM)>0.  (8.2) 


Remark.  By  (8.1 ),  the  bulk  and  interfacial  specific-heats  are  identically 
zero.  Note  also  that,  because  of  (8.2),  the  resulting  interfacial 
free-energy  has  the  form  f(T)  *  e  -  Ts  and  Is  negative  for  large  values  of 
T.  It  is  therefore  tacit  that  -  for  this  model  -  we  are  dropping 
assumptions  (A3)  and  (A6). 

For  convenience,  we  shall  generally  refer  to 

u  *  T  -  Tn  (8.3) 

as  the  temperature.  This  definition  and  (8.1)  allow  us  to  write  the 
interface  condition  (6.19),  In  the  form: 

u*  -hfc/(1-a«,  (8.4) 

where  h  and  a  are  capillarity  lengths  defined  by 

h  -  Tnf(Tn)/L,  a  ■  TMs/L.  (8.5) 


^fs  assumption  is  discussed  by  Mullins  and  Sekerka  ( 1963). 

^Although  wa  make  these  assumptions  for  convenience  only,  we  may  infer  from  remarks  of  Adam 
( 1930],  p.  20  that  for  many  liquids  the  assumptions  concerning  e  and  s  are  reasonable  over  a  wide 
range  of  temperatures. 
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Thus,  by  (6. 1 8M6.20),  we  base  our  model  on  the  equations: 


AU  *  0, 

q  =  -k,gradu 

in  B,, 

(8.6) 

u«-r*y(i-afc), 

Iqlm  *  (L  -  fee)v 

on  I, 

(8.7) 

vmn  =  0 

on  dl. 

(8.8) 

Remarks. 

(1)  By  (8.2),  (8.3)  and  (8.5),  T  computed  using  (8.7),  has  the  unfortunate 
property: 

T  *  0  for  Lleiki  L/(TMs); 

hence  this  range  of  curvatures  lies  outside  the  range  of  validity  of  our 
model.  Mull  ins  and  Sekerka  1 19631  estimate  that  f(Tn)/L  has  10‘7cm.  as 

order-of -magnitude.  If  we  take  this  value  as  an  estimate  for  e/L,  we  see 
that  this  range  of  curvatures  corresponds  to  radii  of  curvature  <  10"7cm. 
(ii)  It  is  not  difficult  to  verify  (by  simply  reversing  the  steps  of  the 
arguments  given  in  Section  5)  that  solutions  of  the  system  (8.6M8.8)  are 
consistent  with  the  first  two  laws  as  presented  in  Section  5. 

In  view  of  the  preceeding  remark,  under  appropriate  boundary 
conditions  solutions  of  (8.6M8.8)  satisfy  the  growth  relations  (7.3)  and 
(7.4)  as  well  as  the  area  estimate  (7.5).  We  can,  however,  deduce  more 
interesting  estimates.  We  begin  by  writing  the  isothermal 
boundary-condition  (7.2)  in  the  form 

u  *  uQ  on  dB  (u0  -  T0  -  Tm).  (8.9) 

Next, 

vol(B,)*  -  -vol(B2)*,  (8.10) 

and, by (8. 1 ) and (8.3), for  gse or s, 

( JE  )•  =  Lvol(B2)*  *  Tm(  JS  )•,  ( Jg  )•  =  garea(I)* 

B  B  i 

Therefore,  if  we  write 


(so  that  8  >  0),  and  assume  that  p  >  O  jn  the  case  of  an  isothermal 
boundary,  then  (7.3)  and  (7.4),  with  T0  -  Trt  in  the  case  of  an  isolated 

boundary,  yield  the  following 

Growth  conditions. 

(i)  For  an  isolated  boundary, 

vol(B2)*  +  Barea(l)*  =  0,  area(I)*  i  0.  (8. 11 ) 

(11)  For  an  isothermal  boundary, 

u0vo1(B1)*+  uarea(I)*  i  0.  (8.12) 


Appropriate  Initial  conditions  for  the  system  (8.6)-  (8.8)  are: 
B,(0)-Bw,  (8.13) 

with  (B0I,B02)  the  initial  distribution  of  phases.  We  shall  designate  by 

Problem  QS  the  initial-value  problem  defined  by  (8.6M8.8),  (7. 1 )  or 
(8.9),  and  (8. 13). 

Questions. 

1 .  What  is  an  appropriate  weak  formulation  of  Problem  QS? 

2.  Is  Problem  QS  well  posed? 

The  growth  conditions  for  QS  yield  a  uniform  bound  on  interfacial  area, 
and  hence  imply  that  the  interface  cannot  become  too  wild.  However,  as 
we  have  remarked  in  Section  5,  estimates  of  this  nature  do  not,  a-priori, 
rule  out  dendrites  with  null  volume.  To  discuss  this  possibility,  let 

£0  *  vol(B2(0))  ♦  Darea(I(0)),  (8. 1 4) 

so  that  E0  is  essentially  the  initial  energy.  Then  we  have  the  following 

mathematically  trivial,  but  physically  interesting,  consequence  of  the 
growth  conditions. 

Coro  1 1  ary .  For  an  isolated  boundary,  if  the  initial  energy  is  small  enough 
(E0  <  V),  then  dendrites  with  null  volume  cannot  form. 
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Proof.  In  view  of  (8. 1 0)  and  (8. 1 1 ), 


vol(B2)*i  0, 

so  that  for  a  dendrite  with  null  volume  to  form,  voI(B2)-»V.  But  by  (8.1 1 ) 
and  (8. 14), 

vol(B2(t))  i  E0,  (8.15) 

and  the  desired  conclusion  follows.  □ 

Remark.  By  (8. 14),  for  B,(0)  a  sphere  of  radius  r,  the  condition  E0>V 

is  equivalent  to  r  <  313.  Using  the  Mullins  and  Sekerka  [1963]  estimate  of 
I0"7cm.  as  an  order-of -magnitude  for  0  =  e/L,  we  see  that  for  the 
formation  of  dendrites  as  discussed  above  we  must  have  r  <  10_7cm. 
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8.2.  Equilibrium  in  the  quasi-stattc  theory. 

We  now  turn  to  a  discussion,  mostly  formal,  of  the  equilibrium  behavior 
of  solutions  to  Problem  QS.  Consider,  then,  a  solution  of  Problem  QS  and 
suppose  that 

as  t**oo,  an  equilibrium  state  (u^,  B^.B^)  exists.  (8.15) 

Let  lm  denote  the  interface  corresponding  to  (B]-fB2- ).  For  convenience, 
we  consider  separately  the  two  types  of  boundary  conditions. 

A.  Isolated  boundary. 

Let  Eq  i  0  be  the  Initial  energy  (8. 1 4).  Then  by  (8. 1 0)  and  (8. 11 )  it 
seems  reasonable  to  expect  that  (B]w,B2w)  is  a  solution,  or  at  least  a 
local  minimizes  of  the  variational  problem:  minimize 

aread^) 

subject  to  the  constraint 
voKBj.)  ♦  Baread^)  »  E0. 

The  solution  of  this  variational  problem,  W  say,  is  given  in  Theorem  2  of 
the  Appendix;  the  results  are: 

(i)  If  E0*V,  then  W  has  a  solution.  Moreover,  every  solution  will  have 

minimal  interface,  and  hence  will  have  an  interface  of  constant  mean 
curvature  which  meets  dB  orthogonally. 

(11)  If  E0  >  V,  then  W  has  no  solution  in  the  standard  sense.  Minimizing 
sequences  (B^B^y  will  have 

vol(Bln)  0,  area(In)  (E0-  V)/B  >  0.  (8. 1 7) 

Remark.  By  our  convention  (Section  2),  phase  1  is  the  solid  phase.  Thus 
if  the  initial  "energy"  is  large  enough,  then,  for  each  minimizing  sequence, 
the  region  occupied  by  the  solid  phase  approaches  a  “set"  which  has  zero 
volume,  but  non-zero  boundary-area,  a  result  which,  by  (7.6),  seems 
indicative  of  the  formation  of  dendrites.  Note  that  this  phenomenon  does 
not  arise  in  the  general  theory  (cf.  the  theorem  of  Section  3.3).  In  the 
quasi-static  theory  the  bulk  energies  are  constant,  and  thus  sufficiently 
large  energies  can  be  attained  only  with  the  formation  of  a  large 
interface;  this  is  the  reason  for  the  behavior  specified  in  (ii)  above. 


On  the  other  hand,  in  the  general  theory  the  bulk  energies  increase  to 
infinity  with  temperature,  and  thus  large  energies  are  attainable  without 
creating  a  large  interface. 

0.  Isothermal  boundary. 

By  (8. 12)  the  equilibrium  configuration  should 

minimize  UqVoKB^)  ♦  paread^). 

This  variational  problem  has  an  obvious  solution.  For  Uq*  0,  Bw  is 
either  B  or  0.  Further,  since 

vol(Bw)  ♦  vol(B2-)  *  V, 

for  u0*0,  B5-  is  either  B  or  0  according  uQ<0  or  u0>0.  Thus, 

granted  the  validity  of  this  variational  problem,  B  is  ultimately  single 
phase,  the  solid  phase  prevailing  if  u0  <  0,  the  liquid  phase  if  u0  >  0.  As 

would  be  expected,  these  results  are  consistent  with  the  Transition 
Theorem  (Section  3.1). 
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8.3.  Quasf-statfc  theory  with  B  unbounded. 

We  now  allow  B  to  be  all  of  space,  but  assume  that 

one  of  the  regions  B,(t)  Is  bounded  for  all  t  (8. 1 8)1 

(cf.  Section  3.2).  Then  for  B2  bounded  the  constitutive  equations  (8.1 ) 

lead  to  finite  values  for  total  energy  and  entropy.  This  will  be  true  also 
for  Bf  bounded  provided  we  stipulate  that  for  this  case  (8.1)  be  replaced 

by 

V-L.  ^“0, 

S,  -  -L/T„,  S2  -  0. 

In  addition,  we  write  Mxl  and  replace  the  boundary  conditions  by 
either  of  the  following  conditions: 

(I)  body  Isolated  at  infinity: 
q(x,t)-*0  as  r-*-oo, 

(II)  body  Isothermal  at  infinity: 
u(x,t)-»u0  as  r-*». 

Then,  since  u  Is  harmonic  in  a  neighborhood  of  infinity,  we  have  the 
following  estimates  as  r+<*>: 

u(x,t)  -  g(t)  ♦  0(r_1),  q(x,t)  -  0(r"2)  if  Isolated  at  infinity, 

u(x,t) «  u0  ♦  0(r_1),  q(x,t)  -  0(r'2)  if  isothermal  at  infinity. 

These  estimates  allow  us  to  establish  (7.3)  and  (7.4)  under  the  present 
hypotheses:  we  simply  replace  D  in  (5.2)  and  (5.3)  by  a  large  ball  whose 
radius  is  ultimately  allowed  to  tend  to  infinity.  Using  the  extended 
versions  of  (7.3)  and  (7.4),  it  is  a  simple  matter  to  verify  the  following: 


1  Then  the  analysis  or  Sections  5  and  6  extend;  we  need  only  require  that  the  first  two  laws  hold  in 
every  bounded  subregion  of  6.  Note  that  the  contact  condition  (6.6)  is  here  vacuous. 
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Growth  conditions  for  an  unbounded  body.  Let  Bj  be  the  bounded 
phase  region 

(i)  tf  the  body  is  isolated  at  infinity,  then 

(-OVoKBj)**  Barea(I)#  =  0,  area(I)#*0.  (8.19) 

( i  i )  if  the  body  is  isothermal  at  infinity,  then 

H)i+,u0vol(B,)*+  uarea(I)*  *0.  (8.20) 


8.4.  Equilibrium  with  B  unbounded. 

We  assume  that  (8.16)  Is  satisfied  and,  for  convenience,  consider 
separately  the  two  conditions  at  infinity. 

A.  Body  Isolated  at  Infinity. 

For  Bj  the  bounded  region  let 

E0  *  the  initial-value  of  (- 1  VvoKB,)  +  Garea(3B,),  (8.2 1 ) 

so  that,  as  before,  E0  is  essentially  the  Initial  energy.  Here  we  have  used 
the  fact  that 


I  -  3B,. 


For  convenience,  write 


D-B. 


for  the  region  occupied  by  the  bounded  phase  at  equilibrium.  Then  by  (8.19) 
it  seems  reasonable  to  expect  that  D  is  a  solution  of  the  variational 
problem: 


minimize  area(3D) 

subject  to  (- 1  )'vol(D)  +  Oarea(dD)  =  E0. 


(8.22) 


Case  1.  Liquid  phase  bounded  (i  -  2). 

The  solution  is  a  sphere  consistent  with  the  constraint. 

Case  2.  Solid  phase  bounded  (t  *  1). 

In  this  case  (8.22)  has  the  equivalent  formulation 


minimize  vol(D) 

subject  to  -vol(D)  ♦  Garea(dD)  ■  E0. 


(8.23) 


For  E0  ■  0  the  solution  is  D  =  0.  For  E0  <  0  the  solution  is  a  sphere 
consistent  with  the  constraint.  For  E0  >  0  (8.23)  has  no  solution; 
minimizing  sequences  (Dn)  have 


vol(Dn)  0,  area  (3Dn)+  E0/G. 


Again  we  have  an  indication  of  a  dendrite  with  null  volume.  If  we  consider 
an  initially-spherical  solid  of  radius  r,  then,  by  (8.21),  for  E0  >  0  it  is 

necessary  that  r  <  313  (cf.  the  last  remark  in  Section  8. 1 ). 


B.  Body  isothermal  at  infinity. 

Case  1 .  Bounded  phase  superthermal  at  u0. 

In  view  of  the  definition  of  the  term  “superthermal"  (cf.  Section  2.1), 
for  this  case  either  i  *  I  and  u0  >  0,  or  I  ■  2  and  u0  <  0.  For  either 

(8.20)  leads  to  the  conclusion  that,  at  equilibrium,  the  bounded  phase 
should 

minimize  lu0lvol(D)  +  pareaOD); 

hence  D  *  0.  Thus  a  superthermal  bounded  phase  should  ultimately 
disappear. 

Case  2.  Unbounded  phase  superthermal  at  u0. 

Here  the  bounded  phase  should 

minimize  -lu0lvol(D)  +  pareaOD).  (8.24) 

This  problem  is  no  different  than  (3.7)  (provided  we  identify1  b  with 
p/u0),  and  our  conclusions  may  be  inferred  from  (i)  of  the  theorem  of 

Section  3.2.  In  particular,  (8.24)  has  no  solution  (other  than  a  local 
minimum  at  D  =  0);  and  (DnJ  is  a  minimizing  sequence  if  and  only  if  it  is 

consistent  with  the  area-volume  limit  and  has  corresponding  volumes 
which  tend  to  infinity.  Examples  of  minimizing  sequences  exhibiting 
dendritic  behavior  are  given  in  Section  3.2. 


^Note  that,  by  (3.6),  the  constant  b  of  Section  3.2  is  approximately  p/lugl;  in  fact,  for  bulk 
'  free-energies  which  are  linear  in  T.  which  is  the  case  here,  the  relation  is  exact. 

i 

i 


i 


PART  III.  THEORIES  BASED  ON  THE  CAPILLARITY  RELATION. 


9.  Approximations.  The  capillarity  relation. 

By  assumption  (A4),  F,(T)  -  F2(T)  when,  and  only  when,  T  •  TM;  thus, 

appealing  to  (6. 1 7)  and  (2.5):  if  interfacial  energy  and  entropy  are 
neglected,  or  if  the  interface  is  flat,  then 

[E]  =  L  and  T  =  TM  on  I.  (9.1) 

If  interfacial  energy  and  entropy  are  not  negligible,  then  the  interfacial 
temperature  T  will  generally  not  equal  the  transition  temperature,  but 
instead  will  depend  -  through  (6. 15)  -  on  the  jumps  in  energy  and  entropy, 
as  well  as  the  curvature,  energy,  and  entropy  of  the  interface. 

We  now  use  (6.14)  to  derive  a  simple  approximation  for  the  interfacial 
temperature  which  is  valid  when  the  interfacial  energies  are  small.  More 
precisely,  we  assume  that 

f(T)«df0(T),  e(T)  *  de0(T)  (9.2) 

with  d  a  small  parameter.  By  (6. 14), 

F2(T)  -  F,(T)  -fcf(T).  (9.3) 

Clearly,  (9.3)  has  the  form  g(T,d)*0  and  for  d=0  has  the  unique  solution 
T=Tm.  Further,  the  partial  derivative  of  g  with  respect  to  T  at  (TM,0) 

is  S1(Tri)-S2(TM ),  which,  by  (2.8),  is  nonzero.  Thus,  by  the  implicit 

function  theorem,  (9.3)  has  a  unique  solution  T=T(d)  near  d=0,  and 
T ■Tri+0(d).  We  therefore  write 


U  =  T-T„ 

and  expand  (9.3)  about  d»0;  using  (2.7)  and  the  fact  that  u=0(d),  we  are 
led  to  the  following  estimate,  valid  to  within  terms  of  0(d2): 

u  ~  -  hfc  (9.4)1 

with  h=0(d)  the  capillarity  length2  (8.5),.  We  shall  refer  to  (9.4)  as  the 
capillarity  relation;  this  relation  is  central  to  all  of  what  follows. 

1  Cf .  Delves  1 19741  for  a  derivation  of  (9.4)  based  on  equilibrium  thermodynamics.  In  Delves’s 
derivation,  as  well  as  in  all  others  of  which  I  am  aware,  the  curvature  enters  through  the  classical 
relation  IpHf  with  (pi  the  jump  in  pressure  across  the  interface,  a  quantity  irrelevant  to  the  theory 
discussed  here. 

O 

‘But  here  f  is  arbitrary,  not  necessarily  of  the  form  assumed  in  Section  8. 


Our  next  step  is  to  derive  a  corresponding  approximation  to  the 
interfacial  balance-law  for  energy.  Trivially,  (2.5)  and  (9.4)  yield,  to  0(d), 


IE)  *  L,  (9.5) 

and  combining  (9.5)  with  (6.4)  leads  to  the  following  interface  condition 
for  energy: 

[q]m  *  Lv,  (9.6)1 

again  valid  to  0(d). 

The  remainder  of  the  paper  will  be  concerned  with  models  based  on  the 
approximate  interface  conditions  (9.4)  and  (9.6). 

Remarks. 

(i)  In  the  foregoing  analysis  the  parameter  d  always  appears  in  the 
product  dfe,  and  therefore  it  is  tacit  in  our  scaling  that  dfe  be  small. 
Thus  equations  based  on  the  capillarity  relation  (9.4)  are  probably  not 
valid  for  problems  involving  large  curvatures' 

(ii)  The  approximation  (9.6),  to  terms  of  0(d2)  rather  than  0(d),  is 

[q]  m  »  (L  *  feh(C,  -  C2)  -  fee)v,  (9.7)2 

provided  e  is  constant.  This  relation  with  the  term  fee  omitted  is  common 
in  the  literature.  However,  the  terms  feh(C,  -  C?)  and  fee  are  both  0(d), 

and  to  neglect  one  without  the  other  is  inconsistent.  In  any  event,  since 
L=0(  1 ),  (9.6)  seems  to  be  the  most  rational  approximation  to  use  in 
conjunction  with  (9.4),  except  possibly  when  the  temperature  away  from 
the  interface  (for  example  at  the  boundary)  is  not  close  to  TM. 


’Also,  as  Rogers  (1983)  has  noted,  the  interfacial  temperature  as  given  by  (9.4)  can  have  negative 
values  for  sufficiently  large  curvatures;  In  fact,  for  fc.>t/f(T^).  (Cf.  (i)  of  the  remark  following 

(8.8).) 

^Cf.  Wollkind  [  1979],  who  derives  (9.7)  assuming  constant  specific  heats. 
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10.  Quasi-linear  theory. 

10.1.  Basic  equations.  Growth  conditions. 

We  now  consider  a  model  based  the  interface  conditions  (6.6),  (9.4),  and 
(9.6)  in  conjunction  with  assumptions  of  constant  specific  heats  C,>0  and 

constant  conductivities  ^>0.  Thus,  writing 


U  -  T  -  T 


M  • 


we  are  led  to  the  system 


C,u* »  k,Au, 

q  *  -k,gradu 

in  B,, 

(10.1) 

u--!^ 

Iqjm  -  Lv 

on  I, 

(10.2) 

vm*n  «  0 

on  dl. 

(10.3) 

for  all  time.  As  before,  we  refer  to  u  as  the  temperature.  Also,  in  view 
of  (2.7)  and  (8.5),  It  is  natural  to  suppose  that 

L>0,  h>0.  (10.4) 

Remark.  The  system  (6. 1 8M6.20),  or  the  simpler  version  (8.6M8.8), 
have  the  advantage  of  being  consistent  with  the  first  two  laws.  Therefore, 
for  an  isolated  boundary  these  systems  have  a  natural  conservation  law, 
balance  of  energy,  and  a  natural  Liapunov  function,  the  total  entropy;  and 
similarly  for  an  isolated  boundary.  Because  of  this  we  were  able  to  obtain 
estimates  such  as  the  bound  (7.5)  on  interfacial  area.  On  the  other  hand, 
because  of  the  approximations  involved,  the  system  ( 1 0.  l )-( 1 0.4)  is,  in 
general,  not  consistent  with  the  second  law  (5.3).  Therefore,  results  such 
as  the  global  growth  conditions  (7.3)  and  (7.4)  and  the  area  bound  (7.5) 
must  be  established  anew. 

Our  discussion  of  ( 1 0. 1 H 1 0.4)  is  based  on  the  following  pair  of 
identities: 

-Lvol(B2)*  -  2jC,u*  -  Jqn ,  (10.5) 

B,  aB 

-hLarea(I)*- 2juC,u*  a  LI Igradul I,2  ♦  juqn , 


(10.6) 


where  INI,  Is  the  norm  defined  by 

Lllpll,2  *  kjlpl2  ♦  k2Jlpl2 .  (10.7) 

B,  B2 

To  prove  ( 1 0.5)  and  ( 1 0.6),  note  first  that  (4.3)  with  g  the 
characteristic  function  for  B,,  (4.6),  and  (10.3)  yield 

area(I)*“  Jfev ,  voKB,)*- fv  »  -  vo1(B2)  .  (10.8) 

i  i 

Note  also  that 

Jqgradu  -  -  Lllgradull,2  .  (10.9) 

B 

By  ( 1 0. 1 ),  ( 1 0.2).  (10.8),  (1 0.9),  and  (6. 1 1 ), 

-Lvol(B2)*a  /[q]m  =  Jqn  ♦  2fC,u*, 

1  *  8i  (10.10) 

-hLarea(I)’-  Ju[q]  m  -  -Jq-gradu  ♦  Juqn  ♦  2juCiir, 

IB  3B  Bj 

which  yield  (10.5)  and  (10.6). 

Equations  ( 1 0.5)  and  ( 1 0.6)  simplify  when 

C,  =C2.  (10.11) 

Indeed,  since  u  is  continuous,  (10.5)  and  (10.6)  reduce  to 

(vol(B2)  ♦  C/u}*  -  -  L"1  jqn ,  ( 1 0. 1 2) 

B  dB 

(harea(I)  +  (C/2)Ju2}#  =  -llgradull,2  -L'^uq-n,  (10.13) 

b  ae 

with 

C  *  Cj/L .  (10.14) 

These  identities  and  the -boundary  conditions  (7.1 )  and  (7.2),  with  the 
latter  written  in  the  form  (8.9),  yield,  as  a  consequence,  the  relations 
listed  below;  there  V  is  the  total  volume,  while  um  ■  um(t)  is  the 

mean  temperature: 

V  *  vol(B),  um=  r’ju. 

B 


(10.15) 


Growth  conditions.  Assume  that  Ct  -  C2. 

(i)  For  an  isolated  boundary, 

(voKBj)  ♦  CVum)*  -  0,  (10.16) 

(harea(I)  -  umvol(B2)  +  (C/2)J(u-um)2)*  i  0.  (10.17) 

B 

(it)  For  an  isothermal  boundary, 

(h  area(I)  -  u0vol(B2)  +  (C/2)J(u-u0)2)*  i  0.  ( 1 0. 1 8) 

B 

Proof.  Consider  first  an  Isolated  boundary.  Then  ( 1 0. 1 2)  and  ( 1 0. 1 3) 
reduce  to  (10. 16)  and 

(h  area(I)  ♦  (C/2)Ju2  )•  *  -I  Igradul I,2  .  ( 1 0. 1 9) 

B 

By  (10.15), 

um(Ju)*  -  Jumu •  =J(umu)-  -  um-Ju  *  ( 1  /2)f  (2(umu)*  -  (um2)*},  ( 1 0.20) 

BBS  B  B 

and  If  we  multiply  (10.16)  by  um,  subtract  the  resulting  equation  from 
( 1 0. 1 9),  and  use  ( 1 0.20)  and  the  fact  that  grad  um  ■  0,  we  arrive  at 
(10.17). 

Consider  next  an  isothermal  boundary.  Then  multiplying  (10.12)  by  u( 

and  subtracting  the  result  from  (10.13),  we  conclude,  with  the  aid  of 
(10.1  A),  that  (10. 18)  holds.  □ 

Remark.  In  ( 1 0. 1 7)  and  ( 1 0. 1 8)  the  i  0  may  be  replaced  by 

*  -llgrad(u-um)llj2  i  0.  (10.211 

Appropriate  initial  conditions  for  the  system  ( 1 0. 1 H 1 0.3)  are: 
u(x,0)  *  g0(x)  for  x  in  B,  Bj(0)  =  B0j,  (10.221 

with  g0  the  initial  temperature-distribution  and  (B01,B02)  the  initial 

distribution  of  phases.  We  shall  designate  by  Problem  OL  the 
initial-value  problem  defined  by  the  field  equations  (10.1),  the 
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free-boundary  conditions  ( 1 0.2)  and  ( i  0.3),  the  boundary  conditions  (7. 1 )  or 
(8.9),  and  the  initial  conditions  (10.22). 

Question.1  What  is  an  appropriate  weak  formulation  of  Problem  QL? 

Since  0  i  vol(B2)  i  V,  ( 1 0. 1 6)  yields  a  bound  for  lum(t)l  which  is 
uniform  in  t,  and  this,  with  (10. 17),  yields  bounds  for  area(Kt))  and  for 
the  L2(B)  norm  of  u(t)-u(-,t);  similar  bounds  follow  from  (10.18).  Thus 
the  growth  conditions  have  the  following 


Then  each  solution  of  Problem  OL  has 


Corollary.  Assume  that  C,-C2. 
(on  its  interval  of  existence) 

area(I(t))  and  Hu(t)llL2(B) 


bounded  uniformly  in  t  by  constants  which  depend  only  on  the  data* 

Remark.  The  bound  on  area(I(t))  precludes  the  formation  of  an  interface 
which  is  too  wild.  This  stability  is  a  consequence  of  interfacial 
free-energy,  manifested  in  the  constant  h.  (As  h  tends  to  zero,  the  bound 
on  area(Kt))  tends  to  infinity.) 

The  relations  (10.1 6H  10.18)  can  possibly  be  used  as  a  basis  for  a 
study  of  existence,  uniqueness,  and  stability  for  Problem  QL,  as  they  yield 
a-priori  estimates  for  u  and  for  the  interface.  Such  a  study,  however,  is 
beyond  the  scope  of  this  paper. 


Cr.  Vlslntin  ( 19841,  who  -  for  a  closely  related  problem  -  gives  a  weak  formulation  and  establishes 
global  existence. 

2i.e„  only  on  8,  C.  h,  u0  (for  an  isothermal  boundary),  and  the  initial  data. 
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10.2.  Equilibrium  In  the  quest-linear  theory. 

In  this  section  we  discuss  the  equilibrium  behavior  of  solutions  to 
Problem  QL,  assuming  throughout  that 

q-c* 

Consider,  then,  a  solution  of  QL  and  suppose  that 

as  t-^»,  an  equilibrium  state  (ut-31-lB>)  exists. 

Further,  let  u^  denote  the  corresponding  mean  temperature  and  I*,  the 
corresponding  interface,  and  let 


(10.23) 


E0  -  the  initial-value  of  vol(B2)  ♦  CVum. 

A.  isolated  boundary. 

By  (10.17)  it  seems  reasonable  to  expect  that 

u„  *  constant «  u__. 


and  therefore,  In  view  of  (10.16)  and  (10.17),  the  equilibrium  state  should 
be  consistent  with  the  variational  problem: 


minimize  hCVaread*,)  ♦  vol(B„2  ){voKB-2 )  -  E0). 


(10.24) 


(in  deriving  (10.24),  was  eliminated  using  the  constraint  equation 

obtained  from  (10.16).)  This  problem  Is  nonstandard,  but  easily  solved. 
Since  hCV>0,  If  (10.24)hasam1nimizer  B-2,  then  (B.,  ,B„2)  has 

minimal  interface  In  the  sense  explained  in  the  Appendix,  in  fact,  solving 
(10.24)  is  equivalent  to  finding  a  v  in  [0,V]  that  minimizes 


hCVWv)  +  v(v "  V- 


(10.25) 


In  view  of  the  properties  of  amjn(v)  expressed  In  Theorem  1  of  the 
Appendix,  (10.25)  has  a  solution  v,  and  the  corresponding  solution  of  Hy 
(cf.  the  Appendix)  minimizes  (10.24).  In  particular,  for  E0  >  2V,  B-2  -  B; 
for  E0  <  0,  B-2  «  0.  Thus,  granted  the  validity  of  (10.24),  the  equilibrium 
configurations  will  have  minimal  interface,  and  hence  ^he  interface  will 
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have  constant  mean  curvature  and  meet  dB  orthogonally.  Further,  by 
(10.23),  the  fluid  phase  will  ultimately  disappear  if  the  initial 
mean-temperature  is  low  enough;  the  solid  phase  will  disappear  if  the 
initial  mean- temperature  is  high  enough. 

B.  Isothermal  boundary. 

By  (10.18),  it  seems  reasonable  fo  expect  that  (B1t.,B2-)  will  solve 
the  variational  problem: 

minimize  u0vol(B.. , )  ♦  h  area(  1^). 

This  is  the  problem  discussed  under  (B)  of  Section  8.2,  and  the  conclusions 
are  no  different. 


kZ*.  -V-  *, 


1 1.  Quasi-static  theory. 

11.1.  Basic  equations.  Growth  conditions. 

A  standard  model  for  solidification  is  based  on  the  system 
( 10. 1H  10.3)  with 

C,-C2-0, 

and  hence  on  the  equations 


AU  -  0, 

q  -  -kjgradu 

in  B,, 

(11.1) 

u«-hfc, 

[q)m  =  Lv 

on  I, 

(11.2) 

vm*n  *  0 

on  dl. 

(11.3) 

These  equations,  although  simple,  have  several  interesting  consequences, 
chief  among  them  being  the  instability  of  the  interface  in  certain 
circumstances.1 

We  begin  our  discussion  of  this  model  by  first  noting  the  following 
consequences  of  ( 1 0. 1 6M 1 0. 1 8): 

Growth  conditions. 

(i)  For  an  isolated  boundary  the  Interracial  area  decreases  with  time  while 
the  phase  volumes  remain  constant: 

vo1(B,>  *  0,  area(I)*  i  0.  (1 1.4) 

(ii)  For  an  isothermal  boundary, 

UqVoKB,)*  ♦  harea(I)*  i  0.  (1 1.5) 

Appropriate  initial  conditions  for  the  system  (l  1.1)-  (1 1.3)  are: 

Bf(0)  *  Bw,  (11.6) 

with  (B01,B02)  the  initial  distribution  of  phases.  We  shall  designate  by 

Problem  MS  the  initial-value  problem  defined  by  ( 1 1 . 1 )-( 1 1 .3),  (7. 1 )  or 
(8.9),  and  (11.6). 

Question.  What  is  an  appropriate  weak  formulation  of  Problem  MS? 


^Mullins  and  Sekerka  (1963,1964).  Saaaiso  Wagner  (1956). 

-57- 


As  in  Section  8.3,  it  is  not  difficult  to  extend  our  results  for  bounded 
8  to  situations  in  which  B  is  all  of  space  with  one  of  the  phase  regions 
bounded  for  all  time  (cf.  (8.18)).  In  particular,  we  have  the  following 

Growth  conditions  for  an  unbounded  body.  Let  B|  be  the  bounded 
phase  region 

(i)  For  the  body  isolated  at  infinity ; 


volte)*  *0,  area(I)*iO. 


(11.7) 


(11)  For  the  body  isothermal  at  infinity ; 


(- 1  )Mu0vol(B()*  ♦  h  area(I)*s  0. 


(11.8) 


Remark.  In  (1 1.4)  and  (1 1.5),  respectively,  the  *0  may  be  replaced  by 
=  -h"1 1  Igradul  I,2  SO  and  ■  -llgradull,2  i  0; 

an  analogous  assertion  applies  to  ( 1 1 .7)  and  ( 1 1 .8). 


mm 


1 1.2.  Equilibrium  In  the  quasi-static  theory. 

Consider  a  solution  of  Problem  MS  and  suppose  that 

as  t-**,  an  equilibrium  state  (uOOIB„1,B..2)  exists. 

A.  Isolated  boundary. 

1.  B  bounded. 

In  view  of  (1 1.4),  we  expect  that  the  equilibrium  configuration  should 
have  minimal  Interfacial -area  (as  defined  in  the  Appendix),  and  hence  the 
equilibrium  interface  should  have  constant  mean  curvature  and  meet  dB 
orthogonally. 

2.  B  unbounded. 

Here,  because  of  (1 1.7),  the  bounded  phase  region,  at  equilibrium, 
should  have  minimal  surface-area  when  compared  to  other  regions  of  equal 
volume.  Granted  this:  the  bounded  phase  region  is  spherical  at  equilibrium. 

B.  Isothermal  boundary. 

Since  ( 1 1 .5)  and  ( 1 1 .8)  are  the  same  as  (8. 1 2)  and  (8.20),  the 
variational  problems  describing  equilibrium  states  in  the  current  theory 
are  exactly  the  same  as  those  discussed  in  (B)  of  Section  8.2  and  (B)  of 
Section  8.4,  and  the  conclusions  are  the  same.  In  particular,  the  assertions 
in  (B)  of  Section  (8.4)  concerning  interfacial  instabilities  -  such  as  those 
exhibited  by  the  prickly  ball  -  are  here  fortified  by  the  work  of  Mullins  and 
Sekerka  [1963],  who  established  interfacial  instabilities  for  unbounded  B 
when  the  bounded  phase  is  an  infinitesimal  perturbation  of  a  sphere. 


1 1.3.  Quasi-static  theory  when  one  phase  does  not  conduct  heat. 

Suppose  that  one  of  the  two  phases,  phase  2  say,  is  a  poor  conductor  of 
heat.  If  we  model  this  by  setting  k2  ■  0;  then,  as  is  clear  from  the 

derivation  of  ( 1 1 . 1 )  -  ( 1 1 .3),  the  temperature  in  phase  2  is  indeterminate, 
while  the  temperature  in  phase  1  obeys  the  relations1 


au  =  0  in  B,, 

u  *  -  hfc,  du/dn  =  L,v  on  dB,, 


(11.9) 


where,  for  convenience,  we  have  assumed  that  the  region  B,(t)  occupied 
by  phase  1  is  bounded,  but  that  B  Itself  is  all  of  space.  Here  1,=  L/k,, 
while  du/3n  is  the  outward  normal  derivative  on  dB,. 

Using  steps  analogous  to  those  used  to  derive  (10.5)  and  (10.6),  it  is  not 
difficult  to  verify  that  (1 1.9)  yield  the  growth  conditions: 


voKB,)**  0,  area(dB,)*s  0. 


(11.10) 


Because  of  (1 1.10),  one  might  expect  solutions  of  (1 1.9)  to  stabilize  at 
large  time  with  B,(t)  ultimately  spherical. 

An  appropriate  initial  condition  for  the  system  (1 1.9)  is 

6,(0)  =  B01,  (11.11) 

with  B01  the  region  occupied  by  phase  1  at  t  =  0. 

Remark.  In  this  model  only  one  phase  enters  the  system  ( 1 1.9),  and  the 
notions  of  supercooling  and  superheating  are  extraneous;  for  that  reason 
the  instabilities  discussed  previously  should  not  be  encountered.  This 
conjecture  is  reinforced  by  ( 1 1 . 1 0)  and,  even  more  so,  by  work  of  Duchon 
and  Robert  [  1 984],  who  establish  local  existence  and  uniqueness  for  the 
problem  (1 1.9),  (1  l.l  1)  in  R2  with  Bot  unbounded.2 
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See  also  Caroli.  Caroli,  Roulet,  and  L anger  [  1985],  who  study  the  existence  of  needle  crystals  using 
the  system  ( 1 1 .9)  (in  a  moving  frame  of  reference). 


Appendix.  Some  variational  problems  involving  interfaces. 

Let  R  be  a  sufficiently  nice  subset  of  B,  and  let  denote  the 

corresponding  interface;  that  is,  is  the  portion  of  dR  contained  in  the 
interior  of  B.  Let 

V  =  vol(B), 

and,  for  each  ve[0,V],  consider  the  problem 

*v  minimize  area(lp)  subject  to  vol(R)  ■  v. 

Let  amin(v)  denote  the  corresponding  minimal  area: 

amin(v)  -  inf{  area(lp):  vol(R)  -  v).  ( 1 ) 

Theorem  1 . 1  For  each  we  [0,  VI,  Hv  has  a  solution ;  and  for  each  solution 

the  corresponding  interface  has  constant  mean  curvature  and  meets  6 >B 
orthogonally.  Moreover,  amjn(v)  is  continuous  on  [0,V],  strictly  positive 

on  (0,V),  and  zero  at  v  =  0,V. 

For  convenience,  we  say  that  a  partition  (R,B\R)  has  minimal 
interface  provided  R  solves  Hv  for  v  *  vol(R). 

Let  J3  >  0  and  q  i  0  be  given.  Consider  the  problem 

Wq:  minimize  areaCl^  subject  to  vol(R)  ♦  Baread^)  *  q. 

Theorem  2.2 

(i)  For  q  i  V,  Wq  has  a  solution,  and  each  solution  R  solves  Hv  for 
some  v*lO,Vl. 

(ii )  For  q  >  V,  Uq  does  not  have  a  solution  in  the  standard  sense,  if  (Rn) 
is  a  minimizing  sequence  with  {In}  the  corresponding  sequence  of 
interfaces,  then 

vol(Rn)  -+  V,  13areaUn)  ■+  q  -  V  >  0, 

so  that  B\Rn  approaches  a  set  with  zero  volume,  but  nonzero  surface 
area. 

'Cf..  Massari  and  Pape  1 1974);  Glusti  [  1981 );  Gonzalez,  Massari,  and  Tamanini  [  1983];  Gurtin  ( 1985). 
^This  solution  of  Problem  Wq  is  due  to  F.  Almgren  (private  communication). 


Proof.  We  can  clearly  write  Wq  in  the  form 

maximize  voKR)  subject  to  vol(R)  +  BareaOlp)  =  q. 

Thus,  since 

amln(v°l(R))  i  areaClp), 
this  motivates  the  problem 

maximize  v  subject  to  v  +  amin(v)  i  q,  0  i  v  i  V.  (2) 

This  problem  clearly  has  a  solution.  For  q  i  V,  the  solution  v  has 
v  ♦  amln(v)  *  q,  and  it  is  not  difficult  to  to  verify  that  any  solution  R  of 

Hv  solves  Wq.  This  yields  assertion  (i)  of  the  theorem.  For  q  >  V,  the 
solution  is  v  =  V,  and  this  yields  assertion  (ii).  □ 


Appendix.  Scaling  the  quasl-1  inear  theory.  Discussion. 

Insight  Into  the  approximations  underlying  the  theories  developed  in 
Sections  10  and  11  can  be  obtained  with  the  aid  of  appropriate  scalings. 
Consider  the  length  and  time  scales  defined  by 

l-f(Ttt)/L,  t  =  t2L/k,TM, 

and  replace  the  coordinates  x  and  t  by  the  dimensionless  coordinates 
x*  =  x/t,  t*  =  t/t . 

Then,  using  the  dimensionless  temperature,  curvature  and  velocity 
u*  *  (T  -  Tm  )/Tm  ,  It*  *  U,  v*  -  vt/C, 

and  the  dimensionless  parameters 

d*  =  TMC,/L,  k*-k2/k,,  c*  =  C2/C,, 

we  can  render  ( 1 0. 1 M 1 0.3)  dimensionless;  the  result,  with  the 
superscript*  omitted,  is 


du  -  AU 

inB,, 

dcu  *  kAu 

in  B2, 

(1) 

u  =  -fc, 

(gradu, 

-  k  gradu2)-m  =  v 

on  I, 

(2) 

vmn  *  0 

on  dl. 

(3) 

Here  time  derivatives  are  with  respect  to  t*  spatial  derivatives  with 
respect  to  x*;  Bj  and  I  are  appropriately  scaled,  and  gradUj  denotes 

the  limit  of  gradu  as  I  is  approached  from  Bj. 

Remarks. 

(i)  The  quasi-static  theory  studied  in  Section  1 1  follows,  as  a  formal 
approximation  to  ( 1 H3),  under  the  assumptions: 

d  «  1,  cd/k  «  1,  k  =  0(1). 

Thus  the  assumptions  underlying  the  quasi-static  theory  are  that  the 
specific  heats  Cj  be  small  compared  to  the  ratio  L/T^,  and  that  the 

order-of-magnitude  of  the  conductivities  k,  and  k.2  be  equal 


(ii)  The  theory  of  Section  1 1.3  (for  phase  2  a  nonconductor  of  heat) 
follows  formally  from  ( I  )-(2)  under  the  assumptions: 

d  «  1,  k  «  1. 

(By  (2),  the  second  assumption  renders  the  temperature  in  phase  2 
irrelevant  to  the  determination  of  the  temperature  in  phase  1.)  Thus  this 
theory  is  based  on  two  assumptions:  that  the  conductivity  in  phase  2  be 
small  compared  to  the  conductivity  in  phase  I,  that  the  specific  heats  Cf 
be  small  compared  to  the  ratio  L/Tn. 

(iii)  The  assumptions 

d  =  0(  1 ),  k  «  1 

yield  the  following  quasi-linear  generalization  of  (1 1.9): 
du  =  au  in  Bt, 

u  =  -fc,  m-gradu^v  on  dBr 

(iv)  Generally,  the  length  scale  l  will  be  very  small.  Thus  the  treatment 
with  0  finite  probably  has  physical  relevance  only  when  the  maximum 
diameter  of  0  (before  rescaling)  is  small,  of  a  size  comparable  to  the 
length  scale  L 

(v)  It  might  be  that  problems  with  0  infinite  -  but  not  all  of  R3  -  arc 
important.  Although  we  have  not  considered  this  case,  it  is  clear  that  the 
basic  results  deduced  in  Section  6,  as  well  as  much  of  what  we  have 
established  in  other  sections,  remain  valid  for  such  0. 
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